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1 . INTRODUCTION 


1.1 GPS BACKGROUND 

The Navstar Global Positioning System (GPS) is a satellite-based radio- 
navigation system that will provide extremely accurate position, velocity, and 
time on a worldwide, contlnous basis. GPS is planned to be operational in 
1989, although limited capability will commence by 1987, and test capability 
is available for a few hours each day now. Since GPS is being designed for 
civil as well as military use, the potential number of users is enormous. 

With potential accuracies under 100 meters, GPS has been shown to have appli- 
cations to almost every phase of land, sea, or air operations. In particular, 
GPS may be a major aid for remote area rotorcraft operations that require pre- 
cise navigation capability. Because of its worldwide coverage and consis- 
tently high accuracy, GPS may be the only source of precise navigation and 
guidance in otherwise unlnstrumented locations in which rotorcraft frequently 
operate. 

The operational GPS configuration will consist of 18-21 broadcasting sat- 
ellites (Including 3 active on-orbit spares), a ground control center, and 
several monitor stations. The monitor stations track the signals from all 
satellites as they make their passes, and relay the information to the Master 
Control Station which computes a best fit predicted emphemerls and clock 
model for the next orbit. At the control center, corrections are applied to 
the range data transmitted from the monitor stations to remove deterministic 
biases. These Include ionospheric delay, tropospheric refraction, general and 
special relativistic effects, antenna phase center offsets, earth rotation and 
timetag corrections. The data are smoothed over a 15 minute Interval by edit- 
ing wild measurements and fitting the remaining data to a least squares poly- 
nominal. This set of smoothed measurements is then processed by a Kalman est- 
imator in the control center. The output state vector includes orbital ele- 
ment perturbations, solar pressure estimates, satellite clock bias, drift and 
drift rate, monitor station clock errors, tropospheric residuals, and polar 
wander residuals. In the final state of this process, clock states are propa- 
gated forward and the reference ephemerls is corrected. These become the 
ephemerls and clock predictions of the navigation message for each satellite, 
which are uplinked to the satellites for continuous broadcast. 

The navigation signal is continuously broadcast on a 20 MHz spread spec- 
trum signal centered at the LI frequency, 1575.42 MHz, and the secondary L2 
frequencey, 1227.60 MHz. LI and L2 are the carrier signals for GPS whose fre- 
quency doppler shifts can be measured for velocity determination. For posi- 
tion determination, the carriers have high rate bl-phase shift keyed codes 
superimposed on them, modulating the phase. There are actually two codes on 
the carrier, Coarse/Acquisition (C/A) code with a chip rate of 1,023 MHz, and 
Precise (P) code with a chip rate of 10,23 MHz. The C/A code is short, re- 
peating every millisecond. Each satellite broadcasts a different C/A code 
chosen from a family of 1,023 specified codes which allows for minimum inter- 
ference between C/A signals from the satellites and thus positive satellite 
identification by the user. The P code is a long sequence, repeating every 
280 days, and each satellite is assigned a week-long portion of this se- 
quence. In addition, a low rate (50 Hz) navigation message is modulated on 
the signal which contains the satellite ephemerides, clock modeling parame- 
ters, satellite status, ionospheric propagation delay parameters, complete 
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satellite constellation almanac, and a special message block. The total navi” 
gatlon message comprises 1500 bits. The P code will be highly protected and 
denied from common use via encryption, due to Its military value. Hence, most 
civil users will be limited to use of the C/A code, which Is Intentionally de- 
graded further to assure its lack of value to enemy military operations. This 
degradation, called "Selective Availability", is currently designed to provide 
100 meters, 2-dimenslonal root-mean-square, accuracy. 

User receivers continuously track any four signals from the six to eight 
satellites in view at any time to solve the hyperbolic positioning problem for 
three position coordinates and an unknown user clock bias. User receivers use 
less precise quartz crystal oscillators which accumulate phase and frequency 
offsets relative to the more stable cesium satellite clocks. The user must 
solve for this clock bias, hence the use of a fourth measurement. In addi- 
tion, the doppler shift of the carrier can be measured and used for velocity 
determination. 

User receivers employ delay lock loops and phase lock loops for tracking 
the signals and extracting the coded navigation message to complete the tri- 
angulation-type navigation computations. In addition to residual phase and 
frequency errors from these tracking loops, other sources of ranging errors 
include mls-modeled signal propagation delays, ephemecis errors, multipath 
errors, and the intentional signal degradation imposed by the DoU on the C/A 
code. 


The position location problem is a matter of computing user location 
given the satellite locations and ranges to the satellites. Actually, since 
the range measurements include the unknown user clock bias, the measurements 
are called "pseudoranges". Similarly, range rate measurements (from the dop- 
pler Information) are called "deltaranges" . With four paeudorange measure- 
ments, tlie position solution can be found from the expression: 
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■» unknown user clock bias 


At^ ” other errors manifested by signal delays 

To minimize costs a typical civil set design will employ only one or two re- 
ceiver channels » sequencing through the four available satellite signals to 
gather the required information. The sequential nature of pseudorange mea- 
surements in the receiver lends very well to recursive Kalman filtering, which 
is the usual technique for computing the navigation solution (as opposed to 
attempting direct solution of the above equation). A low-cost C/A code track- 
ing set will typically employ an 8-state linearized (indirect) Kalman fil- 
ter. The filter may use an adaptive fading memory feature to control possible 
filter divergence by Increasing diagonal elements of the state covariance ma- 
trix whenever smoothed pseudorange measurement residuals are large. This also 
permits more accurate solutions when dynamic conditions are very mild, which 
is the usual case during a typical rotorcraft flight. 

For maximum three-dimensional accuracy, the receiver needs to track well- 
spaced satellites on the horizon and overhead. However, since the satellites 
are continuously rising and setting from the user's vantage point, geometry is 
often less than optimum for resolution of the vertical and two horizontal 
axes. Error contributions from non-optimum geometry of the four tracked sat- 
ellites is called geometric dilution of precision (GDOP). GDOP operates es- 
sentially as a factor that degrades ranging accuracy* For example, a typical 
GDOP factor of 4 will degrade a ranging errors of 9 meters to yield a naviga- 
tion error of about 36 meters when the ranges are resolved into the navigation 
coordinates of latitude, longitude, and altitude. 


1.2 DIFFERENTIAL GPS CONCEPT 


Differential GPS is a concept that elltnlnaten some of the common, bias 
errors experienced by conventional GPS. Differential GPS derives its poten- 
tial from the fact that the measurement errors are highly correlated between 
different users (as well as being highly correlated in time, or autocorrela- 
tion). By employing a second GPS receiver with comparison to truth, slowly 
varying, correlated errors can be isolated and eliminated. In addition, de- 
pending on the relative rates, intentional degradation of the C/A signal may 
be eliminated ly differential GPS as well. Measurement errors are also highly 
correlated between satellites for any particular user, but such common errors 
are removed by the conventional GPS solution as they are indistinguishable 
from user clock bias, hence corrupt only that estimate. 

In differential GPS, a receiver is placed on the ground in the area where 
greater accuracy is desired, as illustrated in Figure 1-1. The correlated 
errors chat a receiver experiences (such as satellite ephemeris errors) should 
be common to all users in a relatively close geographical area. If the static 
user can obtain a reliable estimate of his actual error and transmit that to 
the dynamic user, the dynamic user may be able to compensate for a large por- 
tion of his error. ITie ground user can use an accurate, independent survey ov 
even a long-term smoothed static GPS solution as truth from which to differ- 
ence his real-time po'^ition and velocity estimates. 

Figure 1-1 actually illuetrutes one of several possible variations of 
differential GPS, direct telemetry. Direct telemetry involves sending up 
correction information for four or more siitellltea to the on-board receiver. 

A telemetry receiver transfers the corrections to the GPS receiver. A second 
variation of differential GPS is the ground transmitter, or "pseudo-satellite" 
concept, A ground satellite transmitter generates an actual pseudorandom 
noise modulated GPS carrier that is received by a spare or multiplexed channel 
in the on-board receiver. The error correction information, calculated by a 
co-locatod ground-based receiver, is contained in the digitally encoded mes- 
sage. 


A third variation uses on-board translation: the translator concept. On- 
board signal translation is a non- real- time concept that merely translates the 
received frequency to a telemetry frequency which is tracked on the ground, 

The ground station computes the user location from this signal and modifies it 
with Its own corrections from direct signal traedcing. 
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Aside from 'hese variations of differential GPS^ another possible method 
of error reduction Is relative GPS. In this case, there Is no static ground 
receiver and hence no "truth" position and velocity for the correction set. 
This condition would exist if a survey were not available or If the second set 
were not staionary. Even In this case, the errors of the two setti* arc highly 
correlated, and although neither set can Improve Its position/ velocity esti- 
mate with respect to absolute truth, their relative position/ velocity accuracy 
Is greatly improved. 

This report addressee the direct telemetry mode of differential GPS. 
Within this technique, different user correction schemes are feasible which 
will have varying performance, complexity, and operational tradeoffs. The 
simulation developed makes it possible to study these tradeoffs in a controll- 
able environment. 


1.3 NASA DIFFERENTIAL GPS PROGRAM 


The NASA Ames Research Center is conducting a research program to evalu- 
ate differential GPS concepts for civil helicopter navigation. The specific 
flight phase and mission applications of interest are those where other pre- 
cise aids are not available and for which accuracy requirements are demanding, 
beyond that of conventional GPS. Such applications include remote area search 
and rescue, offshore oil platform approach, remote area precision landing, 
police work, and fire fighting. The program includes system design studies, 
piloted simulation to develop procedures and displays, and a helicopter flight 
test program to demonstrate and evaluate a selected concept mechanization [1]. 

1.3.1 Differential GPS 


Differential GPS is a promising solution to meeting the increased accur- 
acy needs of the rotorcraft community in the applications mentioned earlier. 
Besides removing many of the dominent error sources naturally occurring in the 
system, it has the potential for locally removing much of the effect of the 
intentional signal degradation to be imposed by the DOD on the C/A signal for 
national security reasons [2]. Differential GPS used for non-precision 
approach, or, if possible, precision approach, can be independent of ground 
aids in the terminal area. [3]. Thus, it is a reasonably low-cost 
enhancement to terminal area operations, probably being an add-on feature to 
helicopters already configured with conventional GPS. Furthermore, many other 
users have expressed interest in differential GPS (marine, fixed wing, sur- 
vey), so that the basic capability and differential corrector stations may be 
in place for other applications. 

Therefore, NASA is conducting research of the differential GPS technique 
to determine its potential and to scientifically study and tradeoff various 
implementation Issues. Analysis of these Issues requires a firm understanding 
of the underlying signal environment and filtering techniques applicable to 
implementing differential GPS. The selection of general GPS implementation 
alternatives, as well as specific software implementations within a selected 
alternative, are dependent on the application, error characteristics, and re- 
ceiver characteristics. Therefore, it is important to reasonably investigate 
such tradeoffs before committing to successively more expensive flight test 
programs and ultimately system production. NASA’s responsibility to the pub- 
lic is to foster such research so that the systems ultimately offered to the 
rotorcraft community truly do meet user requirements in the most favorable 
manner, as a result of NASA’s thorough study, simulation, and flight testing. 

1.3.2 System Implementation Alternatives 

There are three basic implementation alternatives for differential GPS: 
direct telemetry ("data link") ground transmitter ("pseudo satellite"), and 
translator. The three concepts have varying advantages and disadvantages for 
different applications. Although only the direct telemetry concept was of in- 
terest for the present study, the other alternatives are briefly described in 
this section. 

In the direct telemetry concept, the corrector or reference receiver is 
located at a surveyed site. The corrector receiver solves for a normal GPS 
solution, then differences this solution with the known location and velocity 


(the latter Is zero in the topocentrlc coordinate frame). The residuals, 
called the differential corrections, are telemetered up to the user receiver 
(usually on some non~CPS telemetry frequency). The user receiver then applies 
these corrections as If they perfectly correlate with his own errors in his 
navigation solution. His ultimate solution Is a "differential solution". 
Correction data for four or more satellites are transmitted to the user. 

In the ground transmitter concept, an actual L-band transmitter Is used 
at the corrector site. This transmitter generates an unused (by the satel- 
lite) portion of the GPS signal that Is received by a spare or multiplexed 
channel of the user receiver. The user receiver thereby locks onto and tracks 
the corrector pseudorandom noise code and carrier. As a "navigation message" 
on this signal, the user receives the correction Information. The correction 
data are calculated by the correction system the same as for the direct tele- 
metry concept; a ground-based receiver differences a conventional solution 
with a survey before it is modulated onto the ground transmitter signal. In 
this case, the user also has the option of performing ranging calculations on 
the ground transmitter signal to generate an additional pseudorange measure- 
ment . 


The third concept Involves signal translation to a new frequency and re- 
broadcast. Actually, two possibilities exist with this technique. A non- 
real-tlne version has a translator on the aircraft which merely translates the 
signals to a now frequency and transmits them to the ground control station. 
The ground control station performs the signal acquisition, lock, track, rang- 
ing, and data demodulation of a normal receiver. The common bias distance 
(delay) between the aircraft and ground control station washes out of tlie so- 
lution like any common delay through user clock bias estimation. Then, the 
ground control station computes the user’s location, and either stores it for 
poRt-f light use, or transmits It up to the aircraft. The control station may 
also use it for other calculations, such as glldepath offsets, which are then 
tramimltted up to the aircraft as steering commands. In all of these proces- 
ses, the final solution at the corrector site Is corrected vdth differential 
corrections from direct signal tracking. 

The other version of the translator concept Inverts this process. Here, 
the translator Is on the ground and retransmits all satellite signals to all 
interested users. Tlie users then perform their own traclclng and correction 
data extraction, knov/lng the survey locations of the translator, autono- 
mous on-board receiver provides the normal solution with which the corrections 
are differenced, 

1.3.3 Mavlgation Filter Alternatives 

As stated earlier, the only Implementation studied in this effort was the 
direct telemetry, or data link, concept. Within this technique, there exist 
several alternatives for generating corrector data at the corrector station 
and for Incorporating these data into the user solution by the user's naviga- 
tion filter. 

The ma,-jor tradeoff is between the use of navigation domain correction 
data (latitude, longitude, altitude or earth-centered, earth fixed - ECEF - X, 
Y, Z and velocities) and tiie use of measurement domain correction data 
(psendorangta, deltarange). If navigation domain corrections are used, the 
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corrections must be applied after, or at the output of, the navigation fil- 
ter. If measurement domain corrections are used, they must be applied before, 
or at the Input to, the filter by differencing them with the receiver's output 
raw pseudorange and deltarange measurements. Thus the corrections in the mea- 
surement domain will have the benefit of filtering, not to mention that the 
user filter will be carrying smaller bias errors in the state estimates. 

A second order tradeoff exists in the case of measurement domain trade- 
offs. The simplest method of pseudorange/deltarange correction is to subtract 
the measurement corrections from the pseudorange/deltarange measurement out- 
puts of the receiver. Then the navigation filter state is the diffei «intlal 
corrected measurements . Another way to Implement measurement domain correc- 
tions is to use them as direct measurements into the filter via an extended 
state vector. If both biases and rates are observable and are reasonably 
well-behaved over the cycle times in which measurement corrections are compu- 
ted and incorporated, then estimation should improve their quality. Of 
course, the wisdom of this approach is highly dependent on the nature of the 
corrections and the rate at which they are Incorporated. Obvious tradeoffs 
exist when one considers the Increased size of the extended state vector and 
corresponding computational burden. 

Other tradeoffs are evident in the implementation of the direct telemetry 
concept. The use of adaptive covariance routines for handling dynamics can be 
investigated. The values of the process and measurement noise terms should be 
adjusted to account for different state error and measurement error character- 
istics. The use of internal ionospheric and tropospheric models, which are 
less precise models than equivalent ground-based models, may be eliminated en- 
tirely if the difference between user and corrector-observed delays is small 
(smaller than the differential error model residuals, that is). The value of 
correcting the velocity states is questionable since velocity errors may be 
more uncorrelated than position errors, and are generally higher frequency 
(less autocorrelation). Incorporation of essentially Independent random er- 
rors will degrade the solution, not improve it. 

Several of these tradeoffs are investigated in this report, using the 
simulation DIFFGPS. 
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II. PROBLEM DEFINITION AND SCOPE 


2 . 1 OBJECTIVES 

The objectives of this effort were to develop a differential GPS simula- 
tion facility and to accomplish some preliminary user navigation filter trade- 
offs. The application for the simulation was civil helicopter operations, 
especially remote area or offshore precise navigation operations which would 
benefit from a system such as differential GPS due to the lack of other navi- 
gation aids. 

To support current flight test experiments being conducted by NASA, a ma- 
jor design guideline was to base relative navigation performance assessments 
on the system being tested in flight, the Magnavox Z-Set. The Z~Set is a 
first-generation GPS receiver developed and tested during the Air Force 
Phase I Concept Validation program from 1975 to 1979, The Z-set is a single- 
channel C/A code tracking receiver. Of significance to this simulation devel- 
opment, the navigation processor employs a sequential measurement-processing 
Kalman filter. The filter models three ECEP positions and velocities, user 
clock bias and bias rate terms, and an optional barometric altimeter bias. 
Therefore, this was the basic configuration of the conventional GPS Kalman 
filter modeled in the simulation. 

The overall scope of the simulation development, discussed in Sec- 
tion 2.2, was to provide a sufficiently high fidelity dynamic motion and prop- 
agation environment to enable accurate comparisons of differential GPS imple- 
mentations. The emphasis was on comparison of differential models, not on 
comparison of conventional GPS with the more accurate differential mode. 

Details of the simulation and differential navigation tradeoff objectives 
are discussed in the next two sections, 

2.2 SYSTEM SIMULATION FACILITY 

Since the simulation was to be used for continuing investigations beyond 
the immediate study, particular effort bad to be applied to assuring its flex- 
ibility, ease of operation, and potential for future enhancement of various 
models. Thus a major objective of the simulation development was to enable 
the non-developer operator substantial Insight into and control of the various 
simulation functions. 

A further objective of the simulation development was to provide a "real- 
istic'’ environment with "realistic" navigation performance output. Ultimately 
the output of the simulation would be compared with time-series flight test 
results. Therefore, it was decided to implement the simulation in a Monte 
Carlo fashion. This did not imply an Interest to run numerous repetitions to 
achieve statlatical significance, however. Instead, the operator was given 
extensive control of input random and deterministic processes to allow crea- 
tion of "representative" cases or replication of a particular desired test 
condition. 

Due to the uncertain nature of the Input data and the difficulty of in- 
tuitively predicting the response of various non-linear Kalman filters, output 
analysis capability had to be very flexible in enabling the operator to access 
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all aspects of the problem. The objective of the post-run analysis program 
development was to provide retrieval of all input data including dynavalcs and 
system error models, and to provide access to many output performance parame- 
ters. In this way the analyst could track down the cause/ effect relation- 
ships of various parameters as well as comparisons between filters, models, 
and implementation schemes. 

2.3 NAVIGATION FILTERS 

Although it was not an objective of this effort to precisely emulate a 
conventional Z-Set Kalman filter, the conventional filter was the foundation 
of the differential filter, therefore was an important starting point. Fur- 
thermore, one configuration of differential GPS to be studied was simple, 
post-filter algebraic correction of navigation domain estimates. A conven- 
tional filter was needed to generate those estimates. The conventional filter 
was developed consistent with available documentation of the "Z-Set” implemen- 
tation. 

There being few precedents for differential filters, there was a bit more 
freedom In implementing these navigation filters. The primary objective was 
to compare measurement-domain corrections with navigation-domain corrections, 
as introduced earlier in Section 1.2.3 and described In detail later in Sec- 
tion 4.2. Beyond that, consideration was given to use of pseudorange only 
versus pseudorange and deltarange corrections, the use of extended state esti- 
mation in the user filter to optimally estimate the error states, the effect 
of the distance between the user and corrector on the differential error cor- 
relation, and the effect of the rate at which corrections are incorporated. 
Some of these effects were simulated while others were considered only in ana- 
lytical form. 
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III. DIFFGPS SIMULATION DESCRIPTION 


3.1 SIMULATION METHODOLOGY 

The methodology for the development of the DIFFGPS simulation was de~ 
signed to address the project objectives described In Section 2. Overall, the 
design was driven ly the fact that the simulation would be operated by anal- 
ysts who were not Intimately familiar with Its development or even with the 
various models. This meant that the user Interface, model parameter access, 
and post-run analysis had to be sufficiently high-level and Intuitive, For 
example, the user Interface Involves mostly menu-select steps and echoing of 
the state of various models and parameters. All Inputs are prompted, thus re- 
quiring almost exclusively yes/no responses or menu Index selections. While 
this interface may be occasionally tlrae-conquming for an experienced operator. 
It assures that un-inltlated users can be productive after very short training 
sessions. And, even for the experienced user, this structure prevents impro- 
perly set up runs and initialization errors. 

As stated earlier, the DIFFGPS was designed as a Monte Carlo simulation 
to provide realism of the environment and its non-linear features, and to pro- 
vide direct relevance to the NASA flight test program. However, the operator 
has access to all error model parameters and even the random number seeds to 
enable him to define a "representative" scenario. Many models (such as satel- 
lite and user motion) can be previewed before the simulation run to determine 
if the selected parameters produce the desired results. The post-run analysis 
program allows plotting of error source pseudorange and deltarange error con- 
tributions in addition to navigation output performance parameters, so that 
performance can be related to the input error characteristics that generated 
it. In addition to this extensive access to the environmental factors during 
post-run analysis, any particular Monte Carlo run is readily initiated by re- 
peating the default random number seeds and subsequent time history. 

A definitive design methodology was adhered to down to the module 
level. Satellite and user motion are modeled with sufficient fidelity for the 
order of the Kalman filter state vector. Constellation alternatives are se- 
lectable for up to 24 circular orbit satellites. User motion is produced by 
the operator establishing a "route plan" by specifying either latitude/ 
longitude/altitude or range/bearing/altitude from an initial location. The 
simulation automatically generates accelerating turns, climbs, dives and lin- 
ear speed changes, checking for adequate maneuvering distances in the route 
plan. 


Particular emphasis was placed on the development of Che error models, 
since they explicitly determine navigation performance. Once again, the 
guideline for determining required levels of fidelity was the order of the 
Kalman filter state which estimated and measured rate terms. Care was taken 
to assure that empirically-backed error levels were incorporated. Dynamics of 
the errors due to user-satell.lte and user-corrector geometry was faithfully 
reproduced according to the literative sources to support representative dif- 
ferential GPS performance conclusion. This caused consideration and modeling 
of some unknown relationships, such as spatial and temporal variations of ion- 
ospheric delay. While the magnitudes of the errors may reasonably be argued, 
the relative comparison of different filters should produce very robust con- 
clusions based on the error model fidelity used. This last feature, reliabil- 
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Ity of the comparative performance of the filters » was the major design me- 
thodology driver in developing the various error models. 

The differential GPS navigation algorithms were the primary emphasis of 
the tradeoff study. To enable an efficient and fair comparison of implementa- 
tion techniques under the Monte Carlo driver environment, the simulation was 
designed to run three corrector and five user navigation algorithms simultan- 
eously during any run, although only one corrector was permitted to provide 
the reference solution for a particular run. Therefore, five user receiver 
navigation algorithms would process Identical motion, satellite geometry, and 
error model Inputs during a simulation run. Due to the modular nature of 
DIFFGPS, it was possible to substitute in different sets or variations of fil- 
ters into the five user and two corrector slots prior to the run. Of course, 
it was equally easy to substitute in variations of error models or even satel- 
lite/user motion models If the need arose. Interface specifications were well 
defined for such flexibility. 

The post-run analysis module required the most flexibility in design to 
provide the analyst sufficient observation of relative motion, error model 
outputs, and Kalman filter operation. This was accomplished by providing 
menu-selectable plots, with selectable data periods and scales, for all param- 
eters of interest. The operator can essentially retrieve all significant in- 
put, intermediate and output variables from a large data file produced during 
the run.. A high-performance plotting package and high-resolution graphics 
provide the display medium, including multiple plots per page and multiple 
curves per plot for comparative analysis. 
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3,2 OVERALL DESCRIPTION AND STRUCTURE 


A conceptual llluBtratlon of DIFFGPS is provided In Figure 3~l. Pseudo- 
range and deltarange measurements are formulated from user and satellite mo- 
tion. The error models corrupt these measurements before they are processed 
by the filter. Operator selections determine which differential and 
conventional filters operate during the run. The simulation alternatively 
processes the corrector and user loops, computing motion and error model 
states successively for each filter. 

The first interaction with the simulation involves satellite 
constellation definition. A nominal 6 plane by 3 satellite constellation is 
stored in the routine, and the operator has the option of redefining or 
modifying the constellation. Specification involves selection of Keplerlan 
elements defining circular orbits. The operator can plot elevation and 
a?,lmuth angles as well as dilution of precision terms. 

User motion is defined by specifying waypoints and dynamics parameters. 

A point mass user is assumed, however maneuvers are constrained to typical 
helicopter performance with respect to the allowable accelerations. The oper- 
ator can interact with absolute or relative coordinates of waypoints in defin- 
ing a nominal flight profile. 

Navigation filters nust be selected for each run. The operator selects 
one or more user navigation filters, and one correction filter. The values of 
filter noise terms, Initial valves, differential correction time period, and 
types of differential corrections can be changed from their default values by 
tlie operator. 

Error models can be selectively "switched off" at setup time. Also, ran- 
dom number seeds can be cbariged and time constants and error magnitudes can be 
adjusted. This allows sensitivity analyses of the effects of the various er- 
ror sourcett. 

During the simulation run, the user and satellites are propagated seqen- 
tially to the corrector and user solution epochs for calculation of relative 
ranges and range rates. The selected filters process these vieasurements , out- 
puting navigation solutions. All relevant data are stored in a file for post- 
slmulatJon plotting and analysis. 

After the simulation terminates, the user can use the post-simulation 
analysis routine to access the time-based data and plot various relation- 
ships. Plot selection is menu-driven and provides the normal selections of 
scales, time periods, plot superposition and labeling. 












3.3 SYSTEM SETUP ROUTINES 


3.3.1 Satellite VlBlblllty 

Satellite visibility can be determined before the simulation run. This 
routine allows the user to select a constellation, plot elevation and azimuth 
angles from any selected location to any number of satellites, and compute 
GDOP, PDOP, and HOOP for any four satellites over a selected period of time. 
Since the operator must designate four visible satellites for the entire per- 
iod of simulation, this routine is useful for determining valid sets. 

An example of the constellation definition is given in Figure 3-2. Any 
parameter is editable for the purpose of setting up the desired set of satel- 
lites. From this point the operator can select elevation and azimuth plots. 

Two types of plots are available. The first is elevation vs. time, shown 
in Figure 3-3. Time scales from a few seconds to 24 hours can be selected, 
and a visibility cutoff elevation angle (greater than zero) can be speci- 
fied. Tlie second type of plot is a polar plot of elevation and azimuth angles 
as shown in Figure 3-4. This plot respresents the view by a user looking di- 
rectly up at zenith (hence, at the center of the plot). 

After selecting a candidate set of four satellites, the operator may plot 
GDOP, PDOP, and HDOP over any period of time for the selected satellites. 

This allows the operator to review the geometry of the constellation chosen 
for the simulation run. A plot of a typical GDOP calculation is shown in Fig- 
ure 3-5, 
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Figure 3-2. DIFFGPS Satellite Constellation Definition 
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3.3.2 User Route Planning 


The operator defines a user flight profile before executing the simula- 
tion. A module is available that provides a flexible, interactive route plan- 
ning capability for the operator. 

The user route planning routine allows the operator to select from a file 
of existing plans, modify an existing plan, or create a new plan. A new plan 
is set up by identifying waypoints in either latitude/longitude/altltude or 
range/bearlng/altltude from an initial location. The simulation will automa- 
tically generate lead turns, vertical accelerations, and linear accelerations. 
Default values are used for terms such as bank angle and g's unless the opera- 
tor identifies new values for a particular segment. 

Upon completion of waypoint data input, the legs are checked for dynaiolcs 
validity such as sufficient anticipatory turn distance in a leg, sufficient 
time to reach a new velocity, or sufficient distance to reach a new alti- 
tude. When the user plan passes all checks, it is labeled ‘'Valid” and can be 
executed and/or filed for future use. A flexible editor is available for mod- 
ifying plans or for correcting legs identified as violating certain validity 
checks. An example user route plan is shown in Figure 3-6. The user profile 
can be plotted before execution, as shown by the horizontal and vertical 
profiles in Figure 3-7, 
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Figure 3-6. User Route Plan 
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3.3.3 Navigation Fllte?: Selection and Specification 

Since the simulation Is designed to run up to five user Kalman filters 
and three corrector Kalman filters simultaneously, a flexible filter selection 
and specification routine is included In DIFFGPS. Typical filters would in- 
clude conventional GPS and various versions of differential filters. Correc- 
tor filters would typically be variations of user filters designed to 
accurately estimate the total observed pseudorange and deltarange errors. 

I 

The operator has the option of selecting any combination of filters to 
run, and must select one corrector filter for the differential correction if 
differential user filters are run. Initial state covariance estimates, pro- 
cess noise terms and measurement noise terms for each filter are set to de- 
fault values or can be individually changed by the operator. Figure 3-8 il- 
lustrates the display of Initialization data in the program. 

In addition to filter parameters, the user has the option of selecting 
the adaptive covariance feature of GPS navigation filters. This selection is 
made during filter setup. 
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Figure 3-8 


Navigation Filter Initialization Data 
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3,3.4 Error Model Selection and Specification 

As with Kalman filter setup, the various simulation error source models 
can be selectively activated and adjusted by the operator prior to simulation 
execution. Selection or deselection of the error models provides the operator 
with the ability to test filter response to individual errors or certain 
groups of errors. For example, the operator may choose to "turn off" Selec- 
tive Availability to observe system performance without the dominant effects 
of that error source. 

The operator may also review and adjust certain parameters in the error 
models, such as correlation time constants and the variance of the error mag- 
nitude. Since the simulation is run in a Monte Carlo mode, the operator may 
also review the random number generator seeds and pick new arbitrary values to 
change the time history of a particular error. Figure 3-9 illustrates a pro- 
gram output of error model parameters. 
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Figure 3-9. Simulation Error Model Parameters 
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3.4 SIMULATION RUN ROUTINES 

3.4.1 Satellite Motion 


Satellite motion la Implemented by solution of the orbital equations at 
the required tlmea. Solution times include once for each user filter cycle 
and once for each corrector filter cycle which may be skewed from the user 
cycle. The orbital element values are converted Into satellte ECEF x, y, and 
z values for subsequent differencing with user ECEF position for range calcu- 
lation. Elevation and azimuth angles are computed, and an elevation "visibil- 
ity cutoff" angle Is calculated. 

3.4,2 User Motion 


The user is propagated along the preplanned route at the selected veloci- 
ties. Accelerations, decelerations, turns, climbs and dives are executed ac- 
cording to the stored user route plan and segment data. 

A basic user segment Is constructed from an entry turn subsegment, an 
exit turn subsegment, and a connecting straight subsegment. Turn subsegments 
are one-half of the angle between successive legs In the route plan. Of 
course, segments may be designated with no tuim subsegment. The route plan 
must start with a straight subsegment, and operator-defined dynamics refer to 
end-of-segment turns. 

Accelerations and altitude changes are accomplished on the straight sub- 
segments. Although attainment of turn bank angle Is assumed to be instantan- 
eous, vertical accelerations are defined to establish steady-state climb or 
descent rates. "Tear drop" turns are not Included, but may be defined by a 
Judicious choice of short segments. 

During execution sub-subsegments are defined at each solution epoch. So- 
lution epochs occur at the same times as satellite solution times for the user 
and corrector filters. Subsegment definition through a turn compensates for 
arbitrary turn start/end points and for accurate velocity along the "true" 
turn path. Figure 3-10 shows a plot of the actual horizontal position for an 
executed user flight profile. 

3.4.3 Error Models 


Eight error sources are modeled in DIFFGPS, including satellite computa- 
tion errors, propagation delays, and user receiver-induced errors. These er- 
rors are modeled to a level of fidelity consistent with the order of the state 
vector of the user Kalman filters and consistent with the goal of providing 
realistic analysis of differential GPS. Each model is discussed in detail in 
the following sections. 






3. 4. 3.1 Satellite Ephemerls Errors 


The satellite ephemerls parameters that are carried In the navigation 
message are the result of a large Kalman estimator operated at the Master Con** 
trol Station. Many relevant paramer.ers are estimated from monitor station ob- 
servations, Including the orbital element perturbations, solar pressure con~ 
stants, satellite and monitor station clock drift parameters, tropospheric re- 
siduals, and polar wonder residual states [4]. The nominal state for the 
first-order correction process implemented by the Kalman filter is provided by 
a batch-processed reference ephemerls computed by the Naval Surface Weapons 
Center. 

Major error sources to date In ephemerls determination appear to be solar 
pressure parameters, gravitational constant uncertainty and satellite clock 
errors. However, the effects of these errors on user-observed (pseudorange 
and doppler) errors are complex. This Is because the nominal ephemerls 
prediction, first order correction, and user ephemerls Incorporation 
algorithms are all least-squares or Kalman filter-type processes. User models 
also incorporate harmonic correction terms In the ephemerls and third order 
fits to satellite clock drift characteristics, all of which are piecewise fits 
to the total orbit with recommended usage spanning only one hour. Thus 
residual errors to the complex, system wide computation process are not easy 
to model, although they are undoubtedly highly correlated errors over any 
particular observation period. A further complication in this process Is that 
lack of complete observability In the estimation processes causes satellite 
clock parameter estimation errors and ephemerls estimation errors to be highly 
negatively correlated, thus often canceling each other. Of course, signifi- 
cant error sources common to all satellites (such as gravitational constants 
and solar pressure) will cause a high degree of error correlation between 
satellites, which Is Interpreted by user receivers as a user clock offset and 
therefore will not affect navigation accuracy. 

In spite of apparent difficulty In developing good models for ephemerls 
errors, three conditions provide some reassurance for the simulation task. 
First, the satellite clock and ephemerls drifts are rather well-behaved due to 
the facts that orbital dynamics are well-damped and that the estimation pro- 
cess contains a lot of smoothing. Secondly, the periods of interest In the 
simulation are relatively short with respect to overall orbital periods or 
even data page spans (piecewise fit periods), so that simple models are rele- 
vant. Third, geometry plays a major role In the error observed by the user, 
since some major effects evident in the orbital plane are largely unobservable 
in the user's ranging domain (cross-track and along-track errors have very 
little effect on line-of-slght range error, for example). Thus It was deemed 
reasonable for this simulation to assume that ephemerls errors have small com- 
ponents that are uncorrelated between satellites, and highly deterministic 
(blas-llke) over the time periods of Interest, with modeling concentrating 
primarily on geometric observations. 


The basic ephemerls errors are modeled as random walks In each of the or- 
bital axes, in-track, DT, cross-track, DC, and radial, DR: 

DT(t) - DT(t - At) + N(0,1) 

DG(t) - DC(t - At) + N(0,1) 

DR(t) - DR(t - At) + N(0,1) 

where; 

o_,„, a„_ are the standard deviations of the random walks 

UX DLi UK 

N(0,l) are three independent zero-mean, unity variance Gaussian 

random sequences. 

The values of the c's were selected to be consistent with typical orbital 
drifts in these coordinates for the GPS satellites. Tlieir values reflect the 
observation and estimation process by the monitor stations and control center. 

Next these errors were transformed into values expressed as orbital pa- 
rameter perturbations of orbit radius (nominally circular orbit) , AR , 
anomaly, AA, inclination, Ai, longitude of the ascending mode, AJ2, and an in- 
termediate angle, A6. These relationships are. illustrated in Figure 3-11. 


The expressions for the orbital parameter perturbations are: 


where: 


AR„ = DR 

i? 


AA 


A0 


^ -1 / DT V 

tan ^ Pj^) 


tan ^ ( 


DC 


Rg 4-DH 


) 


Ai = (A6) cos (A + AA) 
Ail = (AO) sin (A + AA) 


Rg “ nominal radius of orbit 


The new total values of the orbital parameters (e.g., A + AA) are then 
used to compute a new satellite-to-user range magnitude. The difference be- 
tween this range and the nominal range is the pseudorange error due to ephe- 
merls errors. 


28 




ORIGINAL PAGE 19 
OF POOR QUALITY 


3. 4. 3. 2 Satellite Clock Errors 

Satellite clocks exhibit typical atomic frequency standard drifts. One 
of the common representations of clock error Is the Allan variance stability 
measure which models flicker noise, white noise, and Integrated white noise 
tenas to describe the overall variance of the clock drift. However, the GPS 
receiver models satellite clock drift, tjj, with a second order polynominal 
function of time since clock correction epoch: 

‘b " 'bo + # ‘"oc> I 


where 

t^^ » clock bias at epoch time 
Af/f *= clock drift at epoch 
AtQg = time since epoch 

D ® clock drift rate at epoch 

The constants In this expression are transmitted In the satellite navlga 
tion message based on the control station's last computation/upload. To the. 
extent that this receiver modeling matches the exact drift of the satellite 
clock, the residual clock drift error will be modified from the nominal Allan 
variance. 


Eurtherraore, common error sources (such as mls-modeled earth gravita- 
tional effects) may introduce common errors Into the clock drift correction 
terms, especially when sorting out epheraerls errors and clock drift effects. 
For example, satellite clock errors and epheraerls error have been observed to 
be highly negatively correlated. Nevertheless, the satellite clock models 
used in the simulation attempted to emulate typical GPS-observed errors oVer 
short period of time. 

Xn the simulation, the satellite clock was modeled by exponentially cor- 
related frequency noise with additional random walk phase noise. The drift 
rate expression was: 






- At) 


-At/T 


SC 




-2 At/T 


SC 


N(0,1) 


The corresponding clock bias expression was: 

e „ (t) + £ p (t - At) 
2 °CB 

where 


£g^(t) = clock delta range error 
T$c = clock time constant 
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0 = clock drift standard deviation 

CiD 

Gg^(t) " clock pseudorange error 

o » clock bias standard deviation 

The satellites had Independent random processes with nominally equal var- 
iances and time constants. Such terms were, however, operator selectable. The 
standard deviations were scaled by a factor of 3.336 to account for conversion 
of nanoseconds to meters. Also, the standard deviation of the deltarange error 
Is used to model the cumulative deltarange error effects of the other error 
sources, and is adjusted accordingly, 

3, 4. 3, 3 Selective Availability 

Selective Availability is the term used to describe the intentional deg- 
radation of the C/A code (and total denial of use of the P code through en- 
cryption) for national security reasons. Although exact methods for implemen- 
ting Selective Availability are classified, the effect is to cause a ranging 
error in the received signal. The nature of the ranging error is controllable 
by DoD, but is currently advertised to produce an error of 100 meters, two 
distance root-mean-square (2 drms). This was a recent (July, 1983) reduction 
from a previously announced level of 500 meters 2 drms. It is conceivable 
that levels may be reduced in the future, and, of course, DoD will retain the 
capability of raising or lowering the accuracy level whenever security needs 
justify a change [2]. 

More complete statistics of the error levels and their effects were re- 
ported recently by Kalafus, et. al., of the DOT Transportation Systems Center 

[5] . Using declassified samples of tracking segments, Kalafus derived the 
mean and variance bias level, and first and second differences of the er- 
rors. The statistics were calculated from 672 90-second segments of pseudo- 
range errors. Unfortunately, the data reported in that study consisted of 500 
meter, 2 drms class accuracy. For the purpose of this study, KUilafus ' statis- 
tical values were adjusted downward to approximate the 100 meter, 2 drms level 
and similar ratios of first and second difference levels. 

The term 2 drms is, by definition, two-dimensional, orthogonal represen- 
tation of llne-of-slght errors. Thus, 2 drms is a function of HDOP and the 
two dimensions. The operational definition of 2 drms comes from the Federal 
Radio Navigation Plan (FRP) which specified the interpretation of 2 drms 

[6] . It is defined as the radius of the circle containing 95% of the navi- 
gation errors. The FRP further states that it can be assumed that circular 
error probable is related to 2 drms by: 

2 drms =2.5 (CEP) 

To calculate the values of pseudorange Selective Availability error, the 
following convention was established. Pseudorange error, e g^, will occur in 
generally non-orthogonal lines-of-positlon (LOP). If the LOPs were ortho- 
gonal, however, and the two axes had standard deviations of o and a , 
then 2 drias would be defined as [7]: * Y 

2 drms = 2 “J a ^ + a ^ ' 

lx y 
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As an example t and Oy may correspond to the standard deviations of the 
latitude and longitude coordinates. 

Actual pseudorange measurement errors are not orthogonal, of course, and 
suffer a geometric dilution of precision in the conversion to orthogonal coor~ 
dlnates. In the horizontal plane only, HDOP Is the relevant factor relating 
llne-of-position errors to latitude/longitude: 

“hORIZ - <“><>'■> »R 


where: 


‘^HORIZ 


a 


R 


standard deviation of the total horizontal position 

I 2 ? 

error, that is, °y . 

standard deviation of the pseudorange error 


Consistant With the definition of HDOP, all pseudorange errors from the 
four satellites are assumed to be equivalent. Using the above expression, the 
2 drms is: 


2 drms = 2 (HDOP) 

K 

Therefore, for an average HDOP of 1.72, a 100 meter 2 drms value would 
yield a pseudorange error of 29 meters. The values used for all Selective 
Availability parameters, based on scaling of the DOT data, are: 


bias: 100 m, 1 o 

rate: .14 in/s, 1 o 

2 

acceleration: .004 ni/s 1 o 

These values suggest, by the low rate and acceleration values, a highly 
time correlated process. Since the presence of acceleration statistics made 
the exact distribution unimportant to the filter performance, a random walk 
acceleration was assumed, with successive parameters derived by first-order 
"integration" of the higher derivatives; 

£ (t) = e (t - At) + 0 - N(0,1) 
s s S 


where: 

N(0,1) = zero mean, unity variance Gaussian random number 

o_ = standard deviation of Selective Availability, = .004 ra/s‘‘ 

O 


e„(t) + e (t - At) 

£ (t) = s At + e (t - At) 

8 / 8 

£ (t) + £ (t - At) 

c (t) s At + e (t - At) 

s I 8 
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Since a user will see some bias level of an on-going Selective Availabil- 
ity error at any particular time, the simulation was started with initial bias 
levels in each of these parameters f The biases were randomly selected from a 
uniform distribution with end points selected to approximate the standard de- 
viations of each parameter: 

ejO) - U(± .004) 

£ (0) = U(± .16) 

E (0) = U(± 165) 

6 

3. 4. 3. 4 Ionospheric Propagation Delay 

The ionosphere, which is an electron layer located roughly between 50 and 
500 kilometers above the earth, causes delays and refraction of radio- 
frequency signals that traverse the layer. Generally, the refraction is 
small, less than .003 degrees for L-band signals above 10 degrees elevation, 
so is Ignored in this simulation [8], The delay is a function of the total 
electron content (TEC) along the signal path. TEC may vary due to several 
factors; time of day, season, geomagnetic latitude, and the state of solar and 
magnetic activity [9]. In addition, small scale irregularities cause scin- 
tillations. These scintillations are more pronounced at the equatorial and 
polar regions and are stronger at low elevation angles. 

For the differential mode, Che spatial correlation of user and corrector 
ionospheric delays was of most interest. Decorrelation in the delay is due 
primarily to the facts that the two receivers will see different elevation 
angles to each satellite, the path traversed will be at different positions in 
the ionospheric layer, and scintillations will be essentially uncorrelated for 
independent paths. The geometry of the differential ionospheric delay problem 
is illustrated in Figure 3-12. 

The following model does not take into consideration the latitude/longi- 
tude of the ionospheric point (Intersection of signal path and ionospheric 
layer). Most receiver ionospheric models calculate the ionospheric point to 
within a degree. The residual error, therefore, is adequately modeled by the 
random terms of the models. 

The basic ionospheric delay is modeled as a function of the local vertical 
path delay plus an additive scintillation term. For the corrector station, 
this error is: 


where: 

= obliquity factor, ratio of slant path length to 
V vertical path length 

Ey(t) = vertical path delay in nanoseconds 

Eg(t) = scintillation delay in nanoseconds 

c = speed of light conversion for nanoseconds to feet 
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The slant path, as derived by Kruczynski, is [8); 

R.R_SIN 

^ “ (R + ~ h) COB ' EL 1-®^" 
e i 

where 

EL^ = elevation of 1th satellite 
Rg “ radius of earth 

Rg « geocentric radius of the bottom of the ionospheric layer 
R,p “ geocentric radius of the top of the ionospheric layer 
h = user altitude 

The vertical path length is simply: 

Ly = Rj - Rg 

The vertical path delay is modeled as an exponentially correlated random 
variable: 

e y(t) = M -J- EV(t) 

where : 


M = constant nominal delay 
EV(t) is the random part; 

- At/t j -2 At/ 1 ' 

EV(t) « EV(t - At) e + Oj V N(0,1) 


where: 

Tj = ionospheric time constant 
Oj = ionospheric standard deviation 

Although M is certainly not constant in reality, for the purposes of a 
finite simulation run time the variations are adequately modeled by the random 
variables of the model. 


The scintillation term is: 
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where 

“ elevation of 1th satellite 
Og » standard deviation of scintillation term 

The above model is used to generate the Ionospheric delay for the correc- 
tor t The user delay is expressed as a function of the corrector delay to ac- 
count for spatial decorrelation (as well as accounting for different elevation 
angles and scintillations). The user delay is; 

where the slant path length tern is a function of user elevation angles. 

For the user: 


-d /D -d^/D 

Ey(t) = M + EV(t) [e ^ ]+EVU(t)[i-e ] 


This equation substitutes an independent random variable, EVU(t), for the 
random variable EV(t), used in the correction ionospheric model. The substi- 
tution occurs with an exponential weighing of the distance, dj, between the 
two receivers. The Independent process is also an exponentially correlated 
random variablf,: 


EVU(t) = EVU(t - At){ 


-At/ X 


lU 


+ a 


lU 




-2At/t. 


lU 


N(0,1) 


where N(0,1) is a zero-mean Gaussian noise sequence, Independent from the 
sequence driving EV(t) for the corrector. The terras and are nominally 
the same as the parameters used in the corrector model. 

The user scintillation term is assumed to be independent from the correc- 
cor scintillation (a reasonable assumption for all but exactly co-located 
users and correctors): 


£gy(t) = N(0,1) 

sin ( EL^^ + .315^) 

where o„ is the same value used for the corrector scintillation. 

The true ionospheric delay error that a receiver sees is the difference 
between the actual delay, as modeled by the above expressions in this simula- 
tion, and the ionospheric model employed in the receiver (see pseudorange mea- 
surement formulation. Section 3.5.4, 1). The differential user will see only 
the residual between his modeling error and the corrector modeling error, as- 
suming both make corrections. Clearly if the two receivers use different 
models, the possibility for increased errors is present. 
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3. 4. 3. 5 Tropospheric Propagation Delay 

The troposphere causes delay and refraction of signals due to the water 
vapor content in the layer. The major portion of the troposphere extends to 
10 kilometers from the surface of the earth, but effects above that level are 
measurable, including part of the stratosphere which extends up to about 
50 kilometers. The majority of signal delay occurs under 30 kilometers. 

Tropospheric delay is a function of the refractivity along the signal 
path, which in turn is a function of local temperature, pressure, and water 
vapor pressure. An expression for local refractivity, N, as a function of 
these parameters is [10]; 

N = 77.6 ^ + 3.73 (lO"^) 
where: 

P = Pressure in millibars 
T = Temperature in degrees, Kelvin 
e Water vapor pressure in millibars 

All three of these parameters, and hence refractivity, are functions of 
altitude. 

Spatial and temporal variations in refractivity are known to exist. Sea- 
son, latitude, local weather, and other factors that affect temperature, pres- 
sure and humidity will affect refractivity. Satellite elevation angle and re- 
ceiver altitude are strong influences on the total delay due to the predomi- 
nence of effects near the surface of the earth [10], However, horizontal 
inhomogeneities are thought to be generally small for reasonably high eleva- 
tion angles, especially at altitudes above one kilometer [11]. 

Since refractivity is the dominant parameter affecting tropospheric de- 
lay, it was decided to use a model of a refractivity- based delay, with a ran- 
dom refractivity expression. The randomness in the refractivity would account 
for local temporal and spatial variations. As with ionospheric delay, the 
correlation between the user's and corrector's observed tropospheric delay was 
of interest. Therefore, the corrector delay was modeled as a function of this 
local refractivity, and the user’s delay was modeled as a function of the cor- 
rector's delay with the refractivity becoming Increasingly Independent as the 
distance between the user and corrector increased. In addition, the magnitude 
of tropospheric delay is reduced from the surface value as a function of 
altitude for the user. 

The basic tropospheric error for the corrector station is modeled as a 
function of the local refractivity and elevation angle: 

4.79 + .00972 N 9-23 

^ - [.00586 - 360)^ + 294] EL/'-" 

T sin EL, sc 1 

c 1 


37 


where: 


ORIGINAL PAGEJS 
OF POOR QUALH^ 


EL^ * elevation of ith satellite 

The constants are empirical fits to data as developed by Altschuler (10), 

Ngg Is the surface refractlvlty at the corrector site and Is modeled as a 
bias term with an exponentially correlated random part: 


N^^(t) « N + REF(t) 
sc SCO 


REF(t) = REF(t - At) e 


- At/ T 


RFT 


•b o, 


RFT 




-2At/x, 


RFT 


N(0,1) 


where: 

N 


SCO 


Nominal surface refractlvlty (Ngj,^ = 313) 


= Refractlvlty time constant 


Orpt “ Refractlvlty standard deviation 
N(0,1) = Zero-mean Gaussian random variable, unity variance 

The user tropospheric delay Is similar to the corrector delay, but with 
refractlvlty expressed as a function of the corrector refractlvlty to account 
for spatial decorrelation (as well as accounting for different elevation, an- 
gles and altitudes). The user surface tropospheric error is: 


T 

s 


4.79 + .00972 N 
su 

sin EL^ 


- {.00586 (N 

su 


360)^ +. 294] 


where: 

+ I^FCt) + DREF(t) {l-e"^^^^] 

o U dU O 

-At/r I -2At/T ' 

DREF(t) = DREFCt - At) e V l~e N(0,1) 


M = N 

‘^SUO ^^sco 

= Refractlvlty time constant 
°DFT ° ^^f’^^ctlvity standard deviation 


38 






N(0,1) *= Zero-mean, unity variance Gaussian random sequence indepen- 
dent from corrector refractivity random sequence 

d “ Horizontal distance between user and corrector 

D *= Spatial decorrelation scale distance 


Tg is the effective tropospheric error. at the surface. Since tropos- 
pheric delay is a function of user altitude, the final user tropospheric error 
must account for this. Again due to Altschuler, the expression for user trop- 
ospheric error is: 



- 1(6.07 (lO"'^) N + .0213} h + ~ 1*58 (lO"^) } h^] 

T e “ ”su 

6 


v/here: 

h « User altitude in thousands of feet 

The resulting tropospheric error that a receiver will see is the differ- 
ence between the error as modeled by the above expressions and the model em- 
ployed in the receiver (See pseudorange measurement formulation, Sec- 
tion 3. 4. 4.1). The differential user will see only the residual between this 
modeling error and the corrector’s modeling error, assuming both employ 
models. Other errors due to the spatial decorrelation will enter as appro- 
priate. If the two receivers use different internal models of tropospheric 
error, the potential for increased differential errors is present. 

3. 4. 3. 6 Multipath 

Signal multipath exhibits mostly noise-like characteristics for rotor- 
craft at altitude, due primarily to small magnitude reflections off the air- 
frame. Vehicle motion usually precludes any steady reflections of substantial 
magnitude. At lower altitudes multipath signals from the ground may be re- 
ceived that are more deterministic in nature, although off most surfaces, and 
again considering vehicle motion, the reflected signals have a slgnificcnt 
amount of variability. For C/A code users, multipath rejection by the re- 
ceiver tracking loops will take place at altitudes above about 1500 feet. 

The multipath is modeled as an exponentially correlated random variable: 


- At/Xj_^ / -2 At/t„' 

Ej^(t) = E^(t - At) e + Oj^ \ 1 - e ” N(0,1) 
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To account for Increased structure in the multipath signal at lower alti- 
tudes » the time constant is modeled as: 

.5 h > 499 meters 

•01(500-h) h < 499 meters 

Although multipath may be highly correlated between satellites, the value 
of in this model is selected to account for the uncorrelated part of the 
multipath signals. 

3. 4. 3. 7 Receiver Noise 

In this simulation, receiver noise is a term applied to the combined ef- 
fects of pseudorange measurement tracking error and quantization effects. The 
value of receiver noise is dependent on several code tracking loop implementa- 
tion characteristics as well as dynamics. For example, non phase-coherent 
correlation adds an error term. Also, code tracking noise is directly propor- 
tional to code tracking loop bandwidths, and received signal strength (or more 
precisely, carrier-to~nolse ratio) affects receiver noise. 

A general expression has been developed for "normalized" mean-square re- 
ceiver noise due to phase jitter, where the noise variance is normalized by 
the code modulation chip width (293.2 meters for the C/A code). This expres- 
sion includes variable "mechanization parameters" which account for the 
variation in noise contribution due to the particular types of tracking loops 
employed [12], The expression ’ for code loop tracking error is; 


where: 
o 

A =■ 

Bn “ 

Bif = 

C/Nq » 

Ki,K2 = 

Phase-coherent delay-lock 
correlation 

Noncoherent delay-lock , 

separate early/late channels 

Noncoherent delay-lock, 
single early/late channel 





(C/NJ 


irras code loop tracking error 
Chip length, ° 293,1 m for C/A code 

Single-sided noise bandwidth of code tracking loop 

Single-sided bandwidth of the correlator bandpass filter 

Carrier signal power to noise ratio 

Implementation parameters for tailoring the expression to the 
particular type of code loop used, for example: 
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In addition to bandwidth and signal strength considerations for receiver 
noise values, signal dynataics (doppler) can seriously affect receiver noise 
II3). By maintaining carrier synchronization in a phase-locked loop, dynamics 
errors in the code loop are eliminated. However, the carrier tracking loop 
will exhibit a phase error when the carrier frequency varies due to relative 
user-satellite motion (frequency doppler shift). This phase error is: 


0 


0] 

n 


where S) ^ frequency rate of change « 2 ti f 


u) natural frequency of carrier loop 


For a second order phase lock loop with damping factor, d; 

.2 


.2 

U) 

n 


(d + l/4d)‘ 


where: 

= single-sided carrier loop bandwidth 
Using a nominal d = .7071, 

•2 -3 D 2 

u) = 3.56 B., 
n L 


The doppler rate of the signal, f , is related to relative acceleration, 

a, by 


f 


a R 
c 


where 

a = relative acceleration 
R = carrier frequency, 1375 MHz 
c = speed of light 

Combining all the above relations, the phase error, which is related to 
the pseudorange error, is 


0 _% 

Jn ~ X 


c 


R a 

c 2 

3.56 
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For both the phase jitter and dynamic phase tracking error terms, nominal 
values of loop bandwldths (single-sided) were assumed as follows: 


Bjp ■ 100 Hz 
■ 20 Hz 

Therefore, the values of receiver noise terms were 

0 " 6.5 m (for C/N >• 34 dB-Hz) 
o 

2 

Ojj *> .0023 a (a in m/sec ) 


For the satellites with a mean velocity of 3874 m/sec and an orbital rad- 
ius of 2.66 (10^) m. the maximum acceleration (radial) observed by a user is' 
less than 1/2 m/sec*^. Therefore, the dominant user-satellite acceleration is 
due to user motion. 


In the simulation, the receiver noise bias term is modeled as an exponen- 
tially correlated zero-mean random variable with a standard deviation, b , 
equal to the value calculated above: 






-At/l 


■f a 




■2Zkt/T 


N(0,1) 


where: 


Epjjg(t) * value of receiver noise bias at time t 
N(0,1) = zero-mean, unity variance Gaussian random variable 


The acceleration-induced term is modeled as a function of velocity dif- 
ference: 


r ..V _ .0023 [R(t) - R(t - At)] 
rnO'' '' at 

where: 

e j-jjp(t) = value of receiver noise dynamics at time t 
R(t) = user-satellite range rate 


Hie total receiver noise term is simply: 

^rn “ ^rnB ^rnD 
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3, 4. 3. 8 User Clock Errors 

The user receivers’ quartz crystal oscillators contribute errors through 
all timing calculations in the receiver. Of course, user clock bias and drift 
rate are terms in the state vector, so it is errors in these estimates that 
contribute to ranging errors. 

The user clock is modeled the same as the satellite clocks, except that 
the terms in the models are given values conslstant with the stability of the 
quartz oscillator. The clock drift rate and bias expressions are; 


“ 4t)e ^**^'^uc + (Jup'^l-e n(0,1) 






At + H(0,1) 


where 


Ey^(t) = clock delta range error 
T^q = clock time constant 
Oyg = clock drift standard deviation 
Gy^(t) => clock pseudorange error 

0 |,g « clock bias standard deviation 

The standard deviations were scaled by a factor of 3.336 to account for 
the conversion of nanoseconds to meters. 


3.4.4 Pseudorange/Peltarange Measurement Formulation 
3.4.4. 1 Pseudorange 

The pseudorange term in the GPS measurement process is a combination of 
true range, user clock bias, propagation errors, satellite errors, and re- 
ceiver errors (for the ith satellite). 

’’“l “ ‘=<'^R1 ■ + "eI + "cl "si + "ll + "t1 "m "n 

where: 

PRj^ = Pseudorange to satellite i ™ 1-4 
c = Speed of light 

tg^ = Time of reception of signal from ith satellite 

tij^ ■= Time of transmission of signal by ith satellite 


I 


b =* True user clock delay 

e = Llne-of-slght component of the difference between satellite 
location and broadcast ephemerls 


e , “ Satellite clock bias at time of transmission 

cl 

= Line*-of“slght component of Selective Availability value 
= Ionospheric delay along the path traversed by the signal 
^Ti ” Tropospheric delay experienced by the signal 
= Multipath-induced error at the receiver 
= Receiver noise and other receiver-induced ranging errors 

Usually a GPS receiver will correct for ionospheric delay with a stored 
model at the time of pseudorange calculation, but in. the simulation this pro- 
cess was moved into the navigation filter where the tropospheric model re- 
sides. 


Since the simulation computed true user and satellite positions at each 
filter epoch, the range magnitude term, 

C(tj^l - tjj_) 


was calculated directly using the appropriate coordinate tranformatlons of 
user and satellite position. All of the other error parameters were outputs 
of the error models discussed in the previous section. 

3, 4. 4. 2 Deltarange 


Deltarange is the terra given to the doppler-derlved pseudorahge differ- 
ence terra calculated by a receiver. Although a rate term could be calculated 
by a first-difference of successive pseudorange terms, the signal carrier fre- 
quency, at 1.6 GHz, provides a much more accurate reference. The implementa- 
tion of deltarange formulation in a GPS receiver is usually accomplished by 
integrating doppier, counting doppler beats over a "deltarange" interval, 
usually a fraction of a pseudorange cycle. Thus the deltarange measurement is 
generally not an instantaneous rate term, but an average over some small 
interval. Furthermore, the deltarange interval does not span the pseudorange 
measurement epoch (the pseudorange epoch is usually at the start or end of a 
deltarange interval), so there is an additional "latency" in the average 
deltarange value calculated. 

Many of the error sources identified in the previous section will cause 
frequency errors in the received signal. However, since frequency and code 
tracking are implemented in essentially independent tracking loops with their 
own peculiar error characteristics, the relative effect of deltarange and 
pseudorange erors on the respective receiver measurement formulation processes 
is not easy to predict. In addition, emplrcal data (and even concurrence on 
analytical results) about error contributions to signal frequency errors dur- 
ing propagation are not readily available. The major agreement is that such 
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effects are "small" [10, 8J. Indeed, many simulations totally ignore such 
effects (8, 14, 6]. 

Considering a lack of good data, it was decided to group moat deltarange 
errors into a single random term in receiver' noise (although receiver noise 
also Included a deterministic pseudorange term). Where there were obvious 
physical processes that related pseudorange and deltarange errors, such as 
satellite clocks, deltarange errors were modeled explicitly. The model formu 
latlon of deltarange was (for the 1th satellite): 

DR « fl + b + e 
^ £i 

1 - 1-4 


where: 

= Deltarange to ith satellite 

fX “ Average doppler rate, converted to range rate 

b = User clock frequency offset (drift rate) and other random 

range rate errors 

= Satellite clock frequency offset 


Again, since the simulation computed true user and satellite velocity, 
the doppler term was computed directly. However, since, in general, a re- 
ceiver will measure deltarange as an average over the deltarange measurement 
Interval, the time of applicability of the deltarange measurement will be some 
portion of a cycle before the associated pseudorange measurement. In the sim- 
ulation, it is assumed that deltarange measurements are made over one-half of 
the pseudorange cycle, thus the "average" deltarange is accurate approximately 
midway in the deltarange Interval, or one-quarter of a pseudorange cycle 
before the pseudorange measurement epoch. This quarter-cycle latency is 
Implemented by assuming a linear deltarange over the Interval. 

3, 4. 4. 3 Differential Corrections and User Incorporation 

In the differential receiver, one further step beyond simple pseudorange 
and deltarange formulation is necessary. In designs where the differential 
navigation filter process differential measurements, they must be calculated. 

The process begins at the corrector site. If the user is merely to apply 
navigation corrections to his compute latitude, longitude, and altitude, the 
corrector simply computes the difference between a measured navigation solu- 
tion and the survey positlon/statlc velocity. If the user is to compute dif- 
ferences in the pseudorange domain before navigation filtering, then the cor- 
rector must compute pseudorange/deltarange corrections prior to transmitting 
the differences. Of course, in this case the survey location must be expres- 
sed as line-of-slght ranges to the satellites of Interest. The computation of 
corrector station errors is described in Section ,3.4.6. 


It may be argued that deltarange corrections are of little value* If the 
majority of errors In the true relative velocities are due to receiver errors, 
or due to local external sources, or due to short term, spatially Independent 
(over the typical distances between a user and its corrector station source) 
error processes, then the incorporation of differential deltarange errors may 
degrade the Inputs more than helping. Certainly If one error source dominates 
deltarange errors, such as the "rate" of the selective availability error, 
then the use of deltarange corrections can be of value. However, available 
Selective Availability data Indicate that the rate contribution Is small. 
Furthermore, satellite clock drifts are Improving as technology advances, and 
the other error contributions to deltarange are not well understood. There- 
fore, the true value of using deltarange corrections will probably have to 
wait for field test experiments of differential GPS. A simulation at the pre- 
sent time can only have as reliable a conclusion on this subject as the fre- 
quency error modeling is reliable. Nevertheless, both differential pseudo- 
range correction and combined pseudorange/deltarange corrections were computed 
and tested in DIFFGPS. 

The first differential correction employed at the user receiver was the 
navigation correction of ECEF X, Y, Z position and velocity. The output of 
the receiver was corrected by the correction error. For example, In the x co- 
ordinate; 

I'^u = 

where 

= differential ECEF X position 
= user-computed ECEF x position 

Cjj = differential correction to ECEF x coordinate 
computed by the. corrector station 

Although the simulation ensured consistency of tracked satellites, it is an 
important operational consideration that both the user and corrector trade the 
same satellite in this implementation. In the example above, C^ is a function 
of the pseudorange error in each of the four satellites tracked to produce the 
correction term. If the user were tracking a different set of four satel- 
lites, even with only one different such as might occur during a turn with one 
satel.lite "shaded" due to tilting of the user antenna pattern, a completely 
new set of corrections would have to be computed. 

The total possible combinations of four satellites from a maximum of five 
to eight in view are; 



(t) “ 

(4] “ 35 
(4) “ 70 
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Calculations show that a large number of combinations may have "accept- 
able” GDOP at any particular time. For example, out of the possible 35 com- 
binations of a 7-ln-vlew constellation, often 15-20 have PDOPs less than 10. 
The possible number of combinations for which corrections may have to be 
broadcast may be for greater than the 5 to 10 necessary far measurement domain 
corrections. 

3.4.5 User Navigation Filters 
3.4.5. 1 Conventional GPS 

The conventional GPS navigation filter was the foundation for all of the 
filters modeled in this effort. The Kalman filter sequentially processes 
pairs of paeudorange/deltarange measurements from each of the four satellites, 
estimating the position, velocity, and user clock phase and frequency error. 
Satellite broadcast location is used in the ^observation and measurement resid- 
ual processes, and models of atmospheric delays provide estimates of these er- 
ror sources. 

The Kalman filter processing is Illustrated in Figure 3-13, The state 
process noise is first calculated as a diagonal matrix, Q, with diagonal ele- 
ments 2 4 

o At 
a 

^11* ^ 22 > ‘133 “ 4 


2 

‘144 “ % 

2 

qSS* ‘166» ‘177 “ °a 
*188 % 


At + 


Ofa 4t 


At 

At 


where; 


0 ^ = acceleration uncertainty 

Oj^ “ user clock bias uncertainty 

o. ■= user clock rate uncertainty 
b 

An optional (operator-selectable) covariance feature is available to pre- 
vent filter divergence during high dynamics, yet allows high gain operation 
during low dynamic maneuvers. The adaptive covariance term is a function of 
measurement projections. It is Implemented by replacing diagonal terms of the 
covariance matrix with adaptively smoothed terms after the covariance propaga- 
tion, according to the Z-Set formulation [15]. The replacement takes place 
as: 


P ~(i, 1) « P“(i, 1) B 
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Figure 3-13. GPS Kalman Filter Processing 
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The term, B, 1b a smoothed summation of the ratio of the measurement re** 
sldual to the measurement variance projection. The variance projection VP, 
depends on the gain weighting factor: 


VP - HPH'^ + R 


The factor, B, is an average of smoothing parameters for both the pseudo- 
range and deltarange: 

B-ffjli [B^O) 

where : 

j “ jth satlllte 

“ smoothing parameter for pseudorange 
®DR " Bmoothlng parameter for deltarange 

The smoothing parameter, Bj^, is: 


BR(t) - .5 - BR(t-At)] SR(t-At) 

The smoothing factor for the deltarange, B^r, is calculated similarly. 

The various parameters used in these expressions are limited to prevent exces- 
sive control by this factor. 

State propagation is performed for use in the measurement observation ma- 
trix computation, since only the error state is propagated in the filter. The 
state is propagated as a velocity dead reckoner over the cycle period of the 
filter. 

Projected measurements are calculated for the purpose of measurement re- 
sidual formulation. The projected pseudorange is calculated from the received 
satellite ephemeris data, Internal models, and state estimates: 


PRN 


K+1 


K+1 


K+1 


+ b + TR +10 


where : 




*K+1 

%+l 

'K+1 

TR 

10 


projected pseudorange to ith satellite at K+1 
measurement time 

1th satellite position at measurement time 
state propagated to measurement time 
clock bias estimate at measurement time 
model-calculated tropospheric delay 
model-calculated ionospheric delay 


The projected deltarange measurement residual is: 


DPRN, 


K+1 


K+1 

:+ 


- 


-4n 


+ Vl 


K+1 


10^) 


.A 
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Error stace update begins with raeasurement residual formulation. For ex- 
ample, the pseudorange measurement residual, PRMR, is: 


PRMRj^^l = - PWIk+1 


where: 


“ pseudorange measurement at K+1 
PRNj^ 4 .j “= projected pseudorange measurement 

The Kalman gain, K, is calculated from the expression; 
\+l “ ^K+1 ^K+1 t^+l\+l^K+l 


where; 


R= fixed raeasurement noise 

Of course, since measurements are processed sequentially, K can be compu- 
ted as a vector. 

The error state and covariance are then updated with: 



V Y 

K+l K+1 



"k«1 ‘■k+1 


Deltarange measurement, update proceeds similarly. Finally, the state is 
rectified for the next; cycle: 


The implementation of this Kalman filter in DIFFGPS follows the Bierman 
formulation [16][17)[18]., This approach uses the upper triangular fac- 
torization of the covariance: 


P = UDU’'^ 

where U is a unit upper triangular matrix and D is a diagonal matrix. 

This was the formulation used in the Magnavox Z-Set [15], 

3. 4. 5. 2 Navigation Differential Correction 

Differential operation employing navigation domain corrections at the 
output state of the filter involves virtually the same Kalman filter as the 
conventional case. After calculation of the rectified state, the x, y, z and 
Jr, & component differential corrections are subtracted from the state. The 
resulting difference is the differential solution for the navigation correc- 
tion filter. 
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The operator also has the option in this case of applying only position 
corrections. The velocity corrections are set to zero and only the position 
states are corrected. 

3. 4. 3. 3 Measurement Differential Correction 

Measurement correction is the variation of differential GPS that applies 
corrections to the input measurement set of the filter. In this case, the in- 
put to the filter is an "improved" pseudorange/deltarange, with correlated er- 
rors removed from the measurement. 

The filter again operates the same as the conventional filter. Since 
some of the error has been removed from the measurements, theoretically it 
would be possible to "tune" the filter and reduce the measurement noise term, 
R. However, in simulation it was found that the filter was quite robust in 
its operation in either conventional or differential configuration. No clear 
advantage of a different noise vector was observed under the varying error and 
dynamics characteristics. Therefore, to keep the filter similar for the com- 
parison task, it was decided to use identical values for filter parameters in 
all filters. 

3.4.6 Corrector Reference Station Filters 


3. 4. 6.1 Position and Velocity Corrections 

The objective in designing the corrector navigation filter is to have it 
respond, in terms of position and velocity errors, identially to the user nav- 
igation filter to common error sources. The simplest way to implement this in 
the simulation was to use a conventional GPS filter, described in Sec- 
tion 3.4.5. 1, as the corrector navigation filter. 

The navigation solution from the corrector filter included three posi- 
tions and three velocities in the ECEF coordinate frame. The actual corrector 
location was subtracted from the position estimates. The true velocity was, 
of course, zero. The remaining differences were used as corrections for the 
navigation correction type of differential filter. 

The velocity corrections were optional at the selection of the opera- 
tor. If velocity corrections were suppressed, the calculation of errors re- 
mained the same in the corrector filter, but velocity errors were set to zero 
before incorporation in the user filter. 

3. 4. 6. 2 Pseudorange and Deltarange Corrections 

These corrections were used for the measurement correction type of dif- 
ferential filter. In this case, the corrector solution/corrector static posi- 
tion differences were converted from ECEF to line-of-sight pseudorange delta- 
range errors using the corrector's best estimate of the observation matrix, H, 
by; 

z = H X 
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where; ^ 

z = pseudorange/deltarange corrections 
X ” ECEF corrections 

H = observation matrix, direction cosines corrector-to-aatellites 

These corrections are then applied by the differential filter by subtrac- 
ting them from the raw measurements prior to filter processing. As with the 
navigation domain corrections, the operator may select to apply only pseudo- 
range corrections or combined pseudorange and deltarange corrections • 

3.4.7 Simulation Execution 


Simulation processing proceeds through the relevant navigation filter 
cycle periods, moving satellites and the user according to the profiles estab- 
lished during setup. Error models are propagated consistent with the variable 
cycle times, 

A variable in-cycle skew between the user and corrector filter is se.lect- 
able by the operator. If the user filter is processed at times separated 
by At, then the corrector is processed at times At + tsKEW' skew param- 

eter may be selected as any value from 0 to 1. 

In addition to the in-cycle skew, the overall corrector update interval 
is selectable as an integer number of user cycle times. If this option is in- 
voked, differential corrections are updated at the beginning of each corrector 
cycle, with the correction held constant during the cycle. The value applied 
at the beginning of each cycle is the average value of the corrections over 
the corrector cycle. The first correction is applied at an "initial" correc- 
tion time, also selected by the operator. 

The simulation processes the user through segments, subsegments, and sub- 
subsegments. Segments connect waypoints in the profile as selected by the op- 
erator. Subsegments are defined as turns or straights. Within a subsegment 
the user may also execute a pull-up or push-over (vertical acceleration) , a 
climb or dive (vertical velocity), or a linear acceleration or deceleration. 
Sub-subsegments handle sampling of the resultant profile at the user and cor- 
rector solution cycle points. 

During the simulation many parameters are recorded at each filter 
cycle. These include filter outputs, truth data, and several filter parame- 
ters and error model parameters. The file built by this process is accessed 
later in the post-run data analysis and plotting routine. 

3.5 POST-RUN ANALYSIS ROUTINE 

3 . 5 .. 1 Data Analysis and Program Structure 

During the simulation, extensive data are stored at each solution epoch 
over the designated time period of the user route plan. These data are avail- 
able for analysis and plotting after completion of the simulation run. 

The plot menu is shown in Figure 3-14. The major categories include fil- 
ter states before and after measurement incorporation, Intermediate filter 
parameters, error model outputs and position "truth". Mean, standard devla- 
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PLOT MENU t 


REFERENCE - 

USER LAT/LONC PROFILE 1 
USER LATITUDE HISTORY 2 
USER LONCITUOE HISTORY 3 


STATEt COVARIANCE PROPOCATION (ECEFl - 


NAVI-l - X,Y,2, B 

4 

- 7 
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B 
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9 

- 12 

NAVt-1 - TOT - DOT 


13 
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M 
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IB 
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22 
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26 

NAV( + 1 - X,Y,Z, B - DOT 

27 
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31 
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32 

- 35 
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36 
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40 
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DEL(*I • X.Y.Z.B - DOT 

44 

- 47 
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4B 
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32 

' 55 
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36 
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60 
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64 

- 67 
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6S 

' 71 
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72 
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73 
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USER CLOCK BIAS 
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Figure 3-14. DIFFGPS Plot Menu 
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parameters > error model outputs and position "truth”. Mean, standard devia- 
tion and root -mean-square are provided where relevant. The program leaves 
room for expansion in the future. 

3.5.2 Plotting 

Ploting in DIFFGPS accesses the DISPLA plotting package resident on the 
NASA VAX 11/780 computer. This package provides for automatic or manual scp.1- 
ing, label and legend definition, nultlple curves per plot and nultlple plots 
per page. The plots are executed in high resolution graphics on a CRT or can 
be plotted in hardcopy. The post-run analysis routine calculates the correct 
rotations and parameter definitions for executing multiple plots from any sim- 
ulation run. 
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IV. SIMULATION RESULTS 


4.1 CONVENTIONAL GPS 

The conventional GPS solution was computed over the nominal helicopter 
profile described in Sections 3,3.2 and 3.4.2, with all error models, includ- 
ing Selective Availability, in operation. GDOP during this 15-minute period 
was plotted in Figure 3.5 of Section 3.3.1. The GPS receiver performance is 
shown in Figure 4-1. This plot shows the total navigation position error 
which is the difference between the filter best estimate-rectified state and 
truth position. Total error is calculated as the root-sum-square of the three 
axes component errors. The statistics over this period yield a mean error of 

59.1 meters, with a standard deviation of 21.8 meters. The rms error is 
63. C ^meters. 

The major dynamics in the total error are caused by the accelerations of 
the vehicle in the nominal profile: turns, pull-ups, push-overs, and linear 
accelerations/decelerations. To Illustrate these effects. Figure 4-2 presents 
a plot of the latitude of the vehicle on the same time scale as the error plot 
of Figure 4-1. The error dynamics from about 1950 seconds to 2170 seconds are 
seen to be caused by the racetrack turns during this period. The third posi- 
tive "peak" in this series of turns is somewhat higher In magnitude, due to 
the fact that this is the turnout of the racetrack which is executed at a 
1.4 g acceleration as opposed to the racetrack turns which are at about 
1.2 g's. Figure 4-3, a plot of the total velocity error, illustrates the 
relative effects of these turns even better. The other sharp dynamics are due 
to other turns, vertical accelerations and speed changes. 

The slowly varying trend in the total error of Figure 4-1 is due mostly 
to a large bias error in the y axis, as illustrated in Figure 4-4. The mean 
errors in x, y, and z are 14.4 meters, 54.3 meters, and 7,3 meters respec- 
tively. The mean error is due primarily to the Selective Availability errors 
(See Section 4,6.8), which in this satellite configuration and for this par- 
ticular set of Selective Availability errors, happen to manifest themselves in 
the y axis. 

To illustrate the conventional GPS filter performance. Figure 4-5 is a 
composite plot of the four pseudorange measurement residuals. These residuals 
show the filter’s insensitivity to acceleration as errors build rapidly during 
all periods of acceleration. 

As an example of the conventional GPS performance during benign dynamics, 
Figure 4-6 shows the total navigation position error and total velocity er- 
ror. The data are from a two-minute period during the long, straight south- 
west segment in the helicopter profile. A position error of 54.4 meters with 
only a 4.9 meter standard deviation is present. The velocity has a low mean 
of 0.18 meters/second. Again, the major contributor to total position error 
was the y-axis error which had a mean of 52.3 meters. 
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4.2 NAVIGATION CORRECTION, POSITION CORRECTIONS 


The navigation correction differential GPS filter was simulated along 
with the conventional and other filters on the nominal 15~minute helicopter 
profile. Also of significance for Che differential operation was the 
reference corrector filter performance at its static location during the 
simulation. 

The corrector solution Is expected to respond to the various errors simi- 
lar to the airborne solutions, hov/ever without the effects of user dynamics. 
Figure 4-7 demonstrates this, showing the total navigation position error of 
the corrector solution. The error has a mean of 54.7 meters with a standard 
deviation of 11.1 meters. The trend in this curve over the period of the sim- 
ulation should be compared with the corresponding conventional user solution 
in Figure 4-1. Aside from the large error excusions due to user dynamics, it 
can be seen that the long terra trend in the data is Indeed highly correlated 
between the two solutions. 

Figure 4-8 confirms the lack of user dynamics in the corrector solution 
(range rate accelerations due to satellite motion and earth rotation are small 
compared with typical user motion— see Section 3. 4. 3. 7). A quite stable total 
velocity error is present throughout the period with a mean of 0.21 meters/ 
second and standard deviation of 0.15 meters/second. As might be expected, 
the pseudorange measurement residuals are nearly zero-mean and random, with a 
high degree of correlation between the four measurement residuals. A plot of 
these parameters is presented in Figure 4-9. 

The position corrections from this corrector solution were used as the 
navigation corrections for the differential GPS solution as described in Sec- 
tion 3. 4. 4, 3, The navigation position performance Improvement is remarkable, 
as Illustrated in Figure 4-10. Total position error has been reduced from a 
mean of 59.1 meters in the conventional GPS case to 15.3 meters in this 
differential case. This is further illustrated by the x, y, z component 
errors in Figure 4-1 1 . The y coordinate error which dominated the conven- 
tional GPS case errors has been reduced from a mean of 54.3 meters to a mean 
of 2,0 meters. The standard deviation of this error is reduced only from 
22.5 meters to 18.8 meters, hence the less dramatic Improvement in total er- 
ror. Of course, the major contributor to the error variance is the vehicle 
maneuver-caused errors and not steady state biases. The differential correc- 
tions do not improve the maneuver caused errors. 

Velocity errors were only slightly Improved, but, again, these errors are 
dominated by maneuver ing-caused errors and not steady state biases. Fig- 
ure 4-12 shows the velocity errors for the navigation position correction 
case. 
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Figure 4-7. Total Navigation Position Error, Corrector Station 
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Figure 4-8. Total Velocity Error, Corrector Station 
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Figure 4-9. Pseudorange Measurement Residuals, Corrector Station 
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Figure 4—10, Total Navigation Position Error, 
Navigation-Position Correction Differential GPS 
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Figure 4-11. Component Position Errors, 
Navigation-Position Correction Differential GPS 
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Figure 4-12. Total Velocity Errors, Navigation-Positions 
Correction Differential GPS 


4.3 NAVIGATION CORRECTION, POSITION/VELOCITY CORRECTIONS 

This simulation run was identical to the one in Section 4.2 with the ex- 
ception that both position and velocity corrections were applied in the dif- 
ferential solution. The corrector performance was identical to the cases il- 
lustrated in Section 4,2. 

One would expect that the addition of velocity corrections would have 
little effect since the velocity errors are small and relatively uncorrelated 
between different receivers. This is what occured in the simultion, as shown 
in Figure 4-13, Total position errors are essentially unchanged from the pre- 
vious case, and velocity errors actually Increase slightly with the addition 
of velocity corrections. Velocity errors more directly reflect the velocity 
correction errors, and if the velocity corrections are not well correlated 
with the pre-correction user velocity errors, then they are likely to degrade 
rather than improve the post-correction solution. 
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Figure 4-13. Total Position Error and Total Velocity Error 
Navigation-Position/Velocity Correction Differential GPS 
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4.4 MEASUREMENT CORRECTION, PSEUDORANGE CORRECTIONS 

The measurement correction differential filter was simulated along with 
the conventional and navigation correction filters on the nominal IS'inlnute 
helicopter profile. As with the navigation correction filter, the reference 
corrector solution errors are used to Improve the user navigation performance 
(corrector performance is discussed and illustrated in Section 4.2). However, 
in the measurement correction case, the corrector navigation errors are 
converted to equivalent llne-of-slght pseudorange/deltarange errors and 
differenced with the pseudorange/deltarange measurement inputs to the user 
filter. This process is described in Section 3.4.6. 

The total navigation position error is presented in Figure 4~14. The 
mean error is 14.8 meters with a standard deviation of 19,1 meters. This mean 
position error compares with 59,1 meters and 14.3 meters in the cases of con~ 
ventional and navigation correction (position-only corrections) filters, re- 
spectively. The performance Improvement over the conventional solution is 
once again dramatic, and the measurement correction filter performance is 
slightly better than the navigation correction filter. The Improvement is 
probably due to the increased "filtering” afforded by processing the differ- 
enced measurement and corrector information through the Kalman filter. It is 
noteable, however, that the degree of improvement is very modest. 
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Figure 4-14. Total Navigation Position Error, 
Measurement-Pseudorange Correction Differential GPS 
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Figure 4-15 shows the coordinate errors of this case, X, y, and z errors 
in the algebraic correction case of .71 meters, 2.0 meters, and .42 meters, 
respectively, compare with this case of .47 meters, 1.8 meters, and 
.39 meters. As noted in previous comparisons, there Is little change in the 
component error variances. 

The total velocity performance is Illustrated in Figure 4-16. As with 
the navigation correction results, the velocity performance is essentially un- 
changed from the conventional GPS case. 
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4.5 MEASUREMENT CORRECTION, PSEUDORANGE /DELTARANGE CORRECTIONS 

This siinulatlon run was identical to the one in Section 4.4 with the ex~ 
ception that both pseudorange and deltarange corrections were applied in the 
differential solution. The corrector performance was Identical to the cases 
Illustrated in Section 4.2. 

This technique exhibited performance essentially equivalent, although 
slightly inferior to, the pseudorange-only correction filter of Section 4.4. 
This is consistant with earlier observations of the performance comparison of 
position only versus positlon/veloclty corrections in the navigation correc- 
tion filter in previous sections. The total navigation position and velocity 
errors are shown in Figure 4-17, and are very close in magnitude to the 
pseudorange-early correction case. Figure 4-18 shows the x, y, z coordinate 
breakdown of the position errors which are 5Z to lOX worse than the pseudo- 
range-only case. 

To compare measurement correction differential GPS with conventional fil- 
ters under benign conditions, the data during the same two-minute periods as 
used in Section 4.1 were plotted. Figure 4-19 shows the total position and 
velocity errors for the time period 2230 to 2350 seconds. A mean position 
error of 3.2 meters and a mean velocity error of 0.21 meters/second resul- 
ted. Thus the performance improvement during low dynamics, when most large 
errors are due to correlated biases, is more extensive than when dynamics 
errors are significant. 

Figure 4-20 shows the pseudorange measurement residuals for this differ- 
ential filter. The statistics compare roughly with the statistics of the con- 
ventional filter pseudorange measurement residuals, indicating that the con- 
ventional filter is oblivious to the bias errors present in that case as 
should be expected. 
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Figure 4-21, Satellite Ephemeris Errors, User Observation 


76 





ERROR (METERS) ERROR (METERS) ' ERROR (METERS) ERROR (METERS 

0-35 Q.7C1 0.7S 1.2 3.5 1.4 1.6 3,7 1.3 D.7 0.0 G.8 




Of <,iiAur/ 



E^En CR SV 7 D5V)' : 7.3369x10'', 3.0S7i5xL0'", 7.4032X10 ■* 


t^ftBCTCR - FINBPT 

□ - crncn cn 3v to usvy 


V 

1 I 1 1 I I i I 

EPHEN ER SV ID LISV:^ « 3.6263«L0% 3.B03Bm 1O'’, l.eoeeHlD' 


CORPEOTO - PINRFT 
□ - EPHEM CR 5V L2 U5VY 




J L 


J L 


J L 


EPHER ER SV 13 USV)’ J ].3560«L0% 2.tllB«lD^, 1,3953x10' 


Pv- 






OORfe^'TOR - PINRPT 

□ = EPHEn ^5V L7 U5VY - 


I 1 i I I 1 1 ( 1 

ISDO.D 20DD,D 21D0.D 2300.0 Z3DD.D Z4DD.D 2500, D 2600.0 270D.D 

TIME (SEC) 

ERHEO ER 3V 17 U3VJ ; B.8710»<L0‘‘, 3.32S8nlO'®, 6.a72BKlo‘* 


I 

ft ■ 

ft 


Figure 4-22. Satellite Ephemeris Errors, 
Corrector Station Obcervation 
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A. 6 REPRESENTATIVE ERROR ANALYSIS 


In the descriptions of the error models in Section 3. A. 3, the dependence 
of these errors on re] itive motion, receiver location, and random processes 
was Identified. This section shows example plots of simulation outputs of 
each error source to illustrate these effects. 

4.6.1 Ephemeris Errors 

Ephemeris errors are modeled essentially as random walk processes as des- 
cribed in Section 3. 4. 3.1. The conversion of random in-track, cross-track, 
and radial errors to equivalent pseudorange errors allow some difference in 
errors to appear at different locations (see discussion of range decorrela- 
tion, Section 4,7), 

Figures 4-21 and 4-22 Illustrate a single Monte Carlo run of four satel- 
lite ephemeris errors for a user and corrector. The user-corrector separation 
distance varied from about 2 to 20 miles. The initial bias errors are typical 
for the simulation. 

4.6.2 Satellite Clock Errors 


Satellite clocks are modeled by their frequency (rate) and phase (bias) 
offsets as described in Section 3. 4. 3. 2. These drifts are computed in terms 
of pseudorange and deltarange errors. The four clock biases are shown in Fig- 
ure 4-23 for a single Monte Carlo run. Figure 4-24 shows the clock rate for 
the clock from satellite number 12 during that same run. The fact that all 
biases start from zero is inconsequential since each filter estimates and 
tracks the clock bias, 

4.6.3 Ionospheric Propagation Errors 

The ionospheric propagation delay model, discussed in Section 3, 4. 3. 4, is 
a complex function of user-satellite orientation with respect to the ionos- 
pheric shell, elevation angle, user altitude and random effects. The major 
components include a random vertical path delay and a random scintillation 
terra. 


Figure 4-25 shows the nominal value of ionospheric delay as a function of 
elevation angle. The maximum delay is about 95 meters. A typical set of ion- 
ospheric delays for four satellites is shown in Figures 4-26 and 4-27. Sat- 
ellites 10, 17, 12 and 7 are at average elevation angles of 10 degrees, 12 de- 
grees, 25 degrees, and 60 degrees, respectively. User ionospheric delays are 
decorrelated with corrector delays exponentially with separation distance. 

The effects of this are illustrated in Section 4.7 on Range Decorrelation. 
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The eri'or that a user filter will experience is the difference between 
the above ionospheric delays and a user internal, elevation angle-dependent 
model. The residual between the true ionospheric error and the user model 
should be significantly less that the raw ionospheric error. This residual 
was generally quite small, on the order of several meters, for the nominal 
simulation runs. A special case was run with the ionospheric error parameters 
set to a greater degree of randomness (so that user model/true delay residuals 
would be larger) and with range decorrelation distances greatly reduced (so 
that user and corrector stations would observe far more Independent delays). 
The ionospheric errors for the corrector and user from satellite 12 are shown 
in Figure 4-28. Compared with Figure 4-27, it is obvious that the randomness 
of this term has been Increased. 

The effect of these changes on navigation performance was significant. 

The corrector total position' error is shown in Figure 4-29. The mean error 
over this period has been affected moderately, increasing from 54.7 meters 
(see Figure 4-7) to 62.5 meters with the increased ionospheric randomness. 

The standard deviation has more than doubled, from 11.1 meters to 28.5 me- 
ters. Total velocity errors, not shown, also increased significantly for the 
corrector from an rras of 0.26 meters/second to 1.4 meters/second. 

Since the ionospheric errors that caused these changes in the corrector 
navigation error were also selected to be more decorrelated from user observa- 
tions, it is reasonable to expect that the differential solution would de- 
grade. This was the case, as illustrated by the conventional filter and dif- 
ferential filter total position errors in Figure 4-30. The increase in the 
randomness of the conventional solution presented in Figure 4-1 is noteable, 
as is the visually obvious lack of correlation of this plot with the corrector 
solution in Figure 4-29 which was under the same ionospheric conditions. The 
lower plot in Figure 4-30 showns the measurement correction differential fil- 
ter solution, with a 29.2 meter mean error. This yields an error improvement 
of 53% versus 75% for the previous result using more benign, correlated ionos- 
pheric errors. The results reflect the inability of the receiver-based ionos- 
pheric model to reflect the randomness of the ionospheric delay as well as the 
inability of the differential technique to correct for less correlated errors. 
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Figure 4-28, Ionospheric Delay for User and Corrector, 
Satellite 12, Increased Randomness 
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Figure 4-29. Total Navigation Position Error, 
Corrector Station, Increased Ionospheric Delay Randomness 
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4.6.4 Tropospheric Propagation Errors 

Tropospheric propagation delay is a function of elevation angle, user al- 
titude, and a random refractlvlty term as discussed in Section 3. 4. 3. 5. The 
user refractlvlty is a function of corrector refractlvlty, with their indepen- 
dence Increasing exponentially with distance between them. 

Figure 4-31 shows the nominal tropospheric delay as a function of eleva- 
tion angle. Maximum tropospheric delay at low elevation angles is about 
25 meters. A typical set of tropospheric delays from four satellites during a 
Monte Carlo run of the simulation is shown in Figures 4-32 and 4-33. These 
were plotted from the same case as plotted for the ionospheric delays in the 
previous section. The effects of range decorrelation of tropospheric errors 
can be seen in Section 4.7. 


As with ionospheric delay, the user calculates a nominal tropospheric de- 
lay using an internal model. The error experienced by the filter will be the 
residual between the actual delay and this internal model. For the cases 
shown in this simulation, these residuals were generally on the order of a few 
meters. 
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Figure 4-31. Nominal Tropospheric Delay as a Function of Elevation Angle 
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4.6.5 Multipath Errors 
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Multipath errors were discussed In Section 3. 4. 3. 6. Multipath Is modeled 
as a random process which changes as a function of altitude (at low altitudes 
only). The nominal errors In the simulation runs presented in this report as- 
sumed rather low values of multipath* on the order of 1 meter rms. A typical 
plot Is presented In Figure 4-34. For cases where ground reflections are 
small and user motion results in random aircraft surface reflections* this as- 
sumption should be valid. 

However, other studies have recently indicated that ground reflections 
may be significant for the C/A code [18], To investigate the effects of this 
possibility* the standard deviation of the random sequence for multipath was 
Increased from one meter to fifteen meters. Plots of the revised multipath 
errors for a user and the corrector are shown in Figure 4-35, with resulting 
rms values of about 14 meters In each case. User and corrector multi-ath are 
generated by statistically independent random processes. 

The effect of such Increased multipath error on user navigation perfor- 
mance is shown In the following figures. The simulation run is Identical to 
that reported in Sections 4. 1-4. 5, Including Selective Availability errors. 
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Figure 4-34, Multipath Error* Satellite 12 
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Figure A-36 shows the corrector total position error. This plot can be 
compared with Figure 4.7 which had the smaller value of multipath error. Mean 
errors are reasonably similar, but the standard deviation of the corrector er- 
ror has Increased from 11 meters to 22 meters. Total rms velocity errors In- 
creased from 0.26 meters to 1.5 meters. 

The conventional GPS and measurement correction differential GPS total 
position solutions are shown In Figure 4-37, These plots can be compared with 
the earlier cases presented In Figures 4-1 and 4-17 . Conventional user er- 
rors, dominated by Selective Availability errors, are about the same In magni- 
tude, with a moderate Increase In standard deviation from 22 meters to 29 me- 
ters. However, the differential solution suffers significantly under the In- 
creased multipath, now with a mean error of 33.4 meters as opposed to the pre- 
vlojis error of 14.9 meters. Thus, the overall level of Improvement with dif- 
i'erentlal GPS has gone from 15 % to 46Z with the large Increase in multipath 
error. Note that differential corrections Included uncorrelated multipath er- 
rors when Incorporated Into the user differential filter, thus exacerbating 
the effect of user-only observed multipath. 



Figure 4-36. Total Navigation Position Error, 
Corrector Station, Increased Multipath Magnitude 
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4.6.6 Receiver Noise Errors 


Receiver noise generation was discussed in Section 3. 4. 3. 7, Receiver 
noise is modeled as both pseudorange and deltarange errors. Both are mostly 
random in nature, although the receiver deltarange noise has an acceleration 
dependent term. 

Receiver noise was set to reasonably low levels conslstant with C/A code 
tracking performance. Figure 4-38 showns the receiver noise pseudorange and 
deltarange errors during the nominal simulation run. The nos pseudorange er- 
ror was 6.6 meters while the rms deltarange error was 0,012 meters/second. 
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4.6.7 User Clock Errors 


User clock modeling was discussed In Section 3. 4. 3. 8. The user clock Is 
modeled similar to the satellite clocks, except with the parameters set to 
values representative of a quartz oscillator usually employed in user recelv~ 
ers. Both the phase and frequency drifts of the user clock ate modeled and 
converted to equivalent pseudorange and deltarange errors. 

‘ f 

Figure 4-39 Illustrates the user clock performance for one 15-minute run 
of the Monte Carlo simulation. Like the satellite clocks, the user clock er- 
rors start at zero for any particular run. Since the receiver models the 
clock states, only the perturbation values of the clock errors are signifi- 
cant. 
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4,6.8 Selective Availability Errors 


The modeling of Selective Availability errors is discussed in Sec- 
tion 3. 4, 3. 3. These errors are slowly varying biases appearing in the range 
measurements. An example of a set of four Selective Availability erroro from 
four satellites during a simulation run is given in Figure 4-40, The average 
values of the four errors are -19 metersi -47 meters* -31 meters, and 8 me- 
ters. Obviously such slowly varying parameters would vary greatly from run to 
run in the Monte Carlo simulation. 

Since Selective Availability errors are so dominant among the various er- 
ror sources, it was decided to run the nominal 15 minute helicopter profile 
without the Selective Availability errors to assess their impact on the 
various performance conclusions drawn. 

The performance of the conventional filter without Selective Availability 
is shown in Figure 4-41. The mean total position error is 32,1 meters, com- 
pared with 59,1 meters with Selective Availability in Figure 4-1, Velocity 
performance was virtually unchanged. The component errors, in Figure 4-42, 
show where the improvement came from in the total error. The mean x, y, and z 
errors were -14.8 meters, -18,4 meters, and 5.9 meters, respectively. These 
values compare with +14.4 meters, +54,3 meters, and 7.3 meters in the previous 
case with Selective Availability. The Selective Availability in the x-axis 
had the fortuitous effect of moving the mean error across zero in this case. 
The y axis did the same, but the magnitude was greater. Therefore, although 
the total position error showed a change of about 27 meters, the root-sum- 
square of the mean x, y and z component errors was on the order of 78 me- 
ters. This value is closer to the design statistical 2 drma value of 100 me- 
ters. 
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Like the conventional user, the corrector solution performance Improved 
with the removal of Selective Availability. ' The mean total position error of 
the corrector solution, shown in Figure 4-43, was 26.4 meters, compared to 
54.7 meters with Selective Availability in Figure 4-7. The rms component er- 
ror difference of these two cases was again 78 meters. 

The differential solutions from the navigation and measurement correction 
techniques are shown in Figure 4-44. The relevant cases with Selective Avail- 
ability were shown in Figures 4-13 and 4-17. The performance is nearly iden- 
tical to the previous results, indicating that the differential technique had 
removed all of the degradation of Selective Availability, as might be expected 
with the very stable, bias like quantities modeled. Note that the improvement 
in mean position error with differential GPS is about 53% in the absense of 
Selective Availability, compared with 75% in the earlier case. The percentage 
improvement with Selective Availability is conservative, also, since we saw 
earlier that two of the components had somewhat canceling effect on the exist- 
ing errors. 
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Figure 4-43. Total Navigation Position Error, 
Corrector Station, No Selective Availability 
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Range decorrelation is the term applied to the condition whereby certain 
errors become less correlated between different users as the distance between 
them Increases. The causes of range decorrelation are primarily the 
difference in signal paths and differing user-satellite geometry. The level 
of decorrelation can directly related to the spatial variation of physical 
phenomena such as Ionospheric and tropospheric delay. Unfortunately, very 
little detail is known about such variations. The DIFPGPS simulation was 
designed to incorporate a range-dependent decorrelation effect for all 
relevant error sources, with the range scale factor adjustable. Section 4.6.3 
discussed one such change In the range scale factor for ionospheric delay. 

This section provides a more detailed look at the effect of Increasing 
corrector to user distance on differential GPS performance. The case analyzed 
In this section uses the nominal decorrelation parameters of the simulation. 

The user profile simulated for this analysis was greatly simplified from 
the previous l5-mlnute flight. This flight was at constant 1000 feet alti- 
tude, A short straight segment to the west was followed by a mild turn (30® 
bank angle) and a longer straight to the northwest. The total profile took 
about 1-1/2 minutes at 110 knots. The horizontal profile Is shown in Fig- 
ure 4-46. 



Figure 4-46. Horizontal Flight Profile for Range Decorrelation Study 



The corrector station was located at different locations for each run, 
beginning at the origin of the user profile (at 0 feet altitude), then at 20, 
100, and 500 miles away. The corrector navigation performance la shown In 
Figure 4-47. The total position error had a mean of 50.9 meters. The major 
range and geometry dependent error sources are Ionospheric delay, tropospheric 
delay, and ephemerls error. These are shown for the corrector location to 
Satellite 12 in Figures 4-48, 4-49, and 4-50, respectively. 

Ionospheric, tropospheric, and ephemerls errors for the user In this case 
are shown in Figures 4-51, 4-52, and 4-53, respectively. As might be expec- 
ted, these plots compare quite well with the corrector errors. The total mean 
position error of the conventional filter was 49.5 meters as shown In Fig- 
ure 4-54. The lower plot shows the measurement correction differential posi- 
tion error of 4.9 meters. 

The changes In performance at a corrector-user separation distance of 
20 miles was essentially unchanged as shown In Figure 4-55. Analysis of the 
relevant errors In this case shows that, at this distance, they are still 
highly correlated with the corrector-observed errors. 

At a separation distance of 100 miles, the observed Ionospheric error Is 
only slightly changed, its mean value, increasing about 3-1/2 meters. The 
ephemerls error Is also relatively unaffected. However, the tropospheric er- 
ror is markedly decorrclated from the closer ranges as shown In Figure 4-56. 
The resulting navigation performance for the conventional and differential so- 
lutions is shown in Figure 4-57. The differential solution has started to 
show some degradation at this distance. 

At a corrector-user separation distance of 500 miles, the decorrelation 
of error sources becomes more evident. The Ionospheric error now shows some 
Independence as In Figure 4-58. The tropospheric error, shown In the lower 
part of Figure 4-58, is largely decoupled from the earlier examples. The nav- 
igation solutions reflect this for the conventional and differential total po- 
sition errors, Illustrated in Figure 4-59. The conventional solution still 
has a mean error of 50.8 meters, but the differential solution now Improves 
this only to 9.7 meters. 
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Figure 4-49. Tropospheric Delay, Satellite 12, 
Corrector Station, Range Decorrelation Study 
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Figure 4-52. Tropospheric Delay, Dser at 0 Separation 
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Figure 4-56. Tropospheric Delay Error, User Separation of 100 Miles 
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Table A-l shows the rms total position errors for the conventional and 
measurement correction differential filters at the various corrector-user sep- 
aration distances simulated. Although these results Indicate a minor effect 
of range decorrelation, they also reflect the dominance of the Selective 
Availability errors In these cases which Is' assumed to be mostly unaffected by 
corrector-user separation distance. 


Table 4-1. Effects of Corrector-User Range on rms Total 
Position Error (meters) 


Corrector-User 

Separation 

Distance 


(mi) 


Filter Type 



Conventional 

GPS 

Differential 

GPS 

0 

49.7 

6.9 

20 

49.7 

6.9 

100 

49.7 

7.0 

500 

51.0 

10.5 
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4.8 DIFFERENTIAL CORRECTION UPDATE RATES 

In all of the previous examples of differential solutions, the differen- 
tial corrections were made at the same rate as the user filter cycle, that Is, 
at 1 Hz. Since many of the error eources are slowly varying, It may be Just 
as effective to reduce this update rate, particularly If the data link Infor- 
mation bandwidth Is limited. It was possible, In the DIFFGPS simulation, to 
vary the correction update date by specifying a time of first correction and 
cycle time thereafter. If an update frequency lower them 1 Hz is selected, 
the routine averages the 1 Hz corrector difference outputs between updates, 
using the averaged correction for the differential solution throughout the 
specified Interval. 

Two cases were run for the nominal 15-mlnute helicopter case simulated 
earlier. The first used a 30 second cycle, and the second used a 2 minute cy- 
cle. All other conditions were kept constant. The corrector and conventional 
solutions were unchanged, of course. 

The measurement correction differential solution suffered some degrada- 
tion in these two cases. Figure 4-60 shows the differential solution for the 
30 second update rate case. The maan total position error for this case was 
17.1 meters, which compares with 14.9 meters at the 1 Hz correction rate 
(shown in Figure 4-17). The 2 minute update rate differential solution is 
shown in Figure 4-61. The mean error In this case is 19.3 meters. In this 
latter case, one can see the effect of changing differential corrections. The 
series of four racetrack turns are no longer very consistent, and a larger er- 
ror is observed Just after 2300 seconds before the final turns. The three 
test cases are plotted In Figure 4-62. 

The navigation correction differential filters exhibited the same trends 
as the measurement correction case, and even did better for some lower rate 
update cases, possibly due to the benefits of data smoothing offered the 
updating technique. In general, it is clear that lower update rates have lit- 
tle effect on differential GPS performance for the types of error models em- 
ployed in this simulation. 
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Figure 4-62, Change in Total Navigation Position Error, 
as a Function of Correction Update Rate for Measurement 
Correction-Pseudorange/Deltarange Differential GPS 


4.9 NASA Z-SET FLIGHT TEST PROFILE 


NASA executes a standard racetrack flight profile, with ascents and 
descents, for GPS performance evaluation. Flight test data are available from 
such tests. 

The profile was programmed in the DIFFGPS simulation. The horizontal and 
vertical profiles are shown in Figure 4-63. The path Included an initial 400 
foot-per-mlnute climb, a one-mile radius racetrack turn, acceleration from 60 
knots to 100 knots, a 6 mile straight, a final one-mile radius racetrack turn, 
and a 600 foot-per-mlnute descent. The flight test flight path and simulation 
flight path are shown In Figure 4-64. 

The conventional GPS total navigation error, without Selective 
Availability, is shown in Figure 4-65. The average error was 42.3 meters, 
with an rras of 52 meters. The effects of maneuvering are quite evident at the 
latter stages of the flight. 

Figure 4-66 shows a 24 minute segment, comprising one racetrack, of 
latitude, longitude, and altitude position errors from a NASA flight test of 
the Z-Set. The approximate rms total error from this segment of data was 
55 meters, which agrees very well with the DIFFGPS Simulation results. 

Finally, the racetrack pattern was simulated with the measurement 
correction differential GPS model. The total navigation error from this ,run 
is shown in Figure 4-67. The average error is 27,3 meters, a 35% reduction in 
error from the conventional GPS case. This modest level of performance 
improvement as compared to previously discussed results is due to the 
relatively high percentage of maneuvering time in the racetrack profile. 
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V CONCLUSIONS 


5.1 SIMULATION UTILITY 

The DIFFGPS simulation developed during this study has been shovm to be a 
flexible tool for analyzing various scenarios of CPS^refcrenced civil 
helicopter navigation. The Monte Carlo simulation provides the analyst the 
capability to adjust most aspects of the system, the flight profile, the 
receiver Kalman filter, and the signal propagation environment to assess 
differential GPS performance and sensitivities. Major elements of the 
simulation to which the operator has access Include: 

• Satellite constellation and almanac 

• User route plan and flight dynamics 

• Receiver Kalman filter parameters and differential Implementation 

• System ephemerls and clock error characteristics 

• Signal propagation environment characteristics 

A major development In this simulation was modeling of the spatial and 
temporal variations In tropospheric and Ionospheric delays which are 
significant considerations In differential GPS. 

For output, the simulation Includes a wide choice of plots and statistics 
of system Input, Intermediate, and output variables. This enables the analyst 
to observe GPS error Input /performance output' relationships as well as Kalman 
filter operating characteristics. The specific Implementation of the plotting 
and scaling facilitates visual determination of event disturbances and 
parameter correlation and comparison. 

Beyond the general simulation capability, the initial satellite setup 
routine Is useful In Its own right for evaluation of satellite geometry and 
daily visibility. The constellation can be defined to closely match the 
operational satellite configuration, for example, thus providing a reference 
for flight test and other activities. Similarly, the user route planning 
routine can be used to establish actual flight profiles for test activities. 
Because of these features, the simulation Is an effective tool for use In 
conjunction with actual flight tests for extending test conditions or for 
emulation of field test conditions and potential anomaly Investigation. 

5.2 ERROR MODELS AND EFFECTS 

One of the primary efforts in this project was to accurately model GPS 
error sources for emulation of the differential GPS envlornment. The error 
sources modeled Include: 

• Satellite Ephemerls Error 

• Satellite Clock Error 

• Selective Availability 

• Ionospheric Propagation Delay 

• Tropospheric Propagation Delay 

• Multipath 

• Receiver Noise 

• User Clock Error 


Other mechanization error sources, such as Kalman filter modeling errors, 
were computed explicitly In the receiver navigation filter simulation. 

The dominant error source was, of course, Selective Availability. Total 
errors without Selective Availability were on the order of 30 meters for the 
Z-Set (Including representative dynamics); errors with Selective Availability 
active were typically 100 meters or more. However, Selective Availability was 
modeled as a very low dynamic error and was virtually eliminated from the 
final solution of the differential GPS Implementation. 

Ionospheric and tropospheric errors are of particular Interest to 
differential GPS because of their spatial variation, therefore possibly 
affecting user and corrector differently. Also, since the user receiver 
typically employs Internal models of these errors, there Is a tradeoff, In the 
differential mode, between eliminating these user models and ensuring matched 
models between users and correctors. Nevertheless, Ionospheric and 
tropospheric errors were usually small (a few meters, worst case). 

Tropospheric errors were mostly decorrelated at user-corrector separation 
distances of 100 miles or more, while ionospheric errors did not show much 
independence until at least 500 miles. Differential GPS errors doubled from 
the co-located case, but still saw significant Improvement due to the highly 
correlated Selective Availability errors. 

Multipath is not a significant error contributor at high altitudes since 
the reflection surfaces are a short distance from the GPS antenna (e.g., 
aircraft surfaces) and the dynamics of flight cause them to be mostly specular 
In nature. However, at altitudes below about 1500 feet, the possibility 
exists for ground-reflected or water reflected multipath which may be strong 
and reasonably constant. Since there were no significant test data on GPS C/A 
code multipath, two cases were simulated, one with a l-raeter rms and one with 
about 14 meters rms. The 14 meter rms case caused a 100% increase In 
differential GPS total navigation position error. In the conventional GPS 
case, the effects of Increased multipath were mostly unobservable due to the 
dominance of the Selective Availability errors. 

The other errors were modeled at nominal GPS error budget levels and 
completed the contributions to total system error. Error analysis and 
modeling continues to be a ripe area for investigation in GPS, particularly 
with those processes that have not been observed to any extent in field 
testing. As tests continue, particularly with multiple sets and differential 
tests, these models can be modified to reflect the empirical data. 

5.3 DIFFERENTIAL GPS IMPLEMENTATION 

Differential GPS was implemented In two primary modes in the simulation, 
utilizing navigation domain (unflltered) differential corrections and 
measurement domain (filtered) differential corrections. Additionally, the use 
of pseudorange/posltlon only corrections versus the use of both 
pseudorange/position and deltarange/velocity corrections were studied. 

The major obvious conclusion was that the Improvement afforded by 
differential GPS over conventional GPS was significant no matter what the 
Implementation technique. Over a composite flight profile, the difference in 
differential GPS performance due to the Implementation technique was small. 
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However, the study did not determine the validity of this conclusion for 
special conditions in Isolated applications. Certain remote applications or 
cases in particular phases of flight may identify advantages of one technique 
or another. Within the cases tested and for the filters modeled, differential 
GPS total position error varied by at most 3-'4Z. 

I 

As might be expected, the effects of user-to-corrector reference station 
distance were minimal for distances of 100 miles or less. The range 
correlation scale factors were based on estimates in the literature of the 
spatial variation of relevant parameters. Empirical data on GPS signal 
propagation effects will be valuable in validating estimates of the range 
decorrelation of these effects. 

Differential GPS performance was relatively insensitive to correction 
update rates as well. Increasing the correction period from 1 second to 2 
minutes caused a 30% increase in differential GPS total position error, but 
the overall improvement from the conventional case was still 67%, 

With the models tested, the overall Improvement afforded by differential 
GPS was about 75% over the composite flight profile, and over 90% for a short 
period of no maneuvering. Without Selective Availability, the improvement was 
44% and 87%, respectively. Table 6,1 illustrates these results. 

Table 6.1 Differential GPS Simulation Results 
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5.4 RECOMMENDATIONS FOR FURTHER STUDY 

Throughout the courue of this study, certain areas requiring further 
research became evident. Tliese were: 

• Measurement of data and model Improvement for propagation errors 

• Improvement of filter performance In dynamics (for those applications 
requiring such precision) 

• Investigation of extended-state and other differential filter 
mechanizations 

• Additional filter adaptive features 

• Filter aiding by available on-board sensors 

• Use of tropospherlc/lonospherlc models by the user and corrector 
receivers. 

Investigation of these areas would enhance the fidelity and validity of 
the simulation and elucidate the implementation tradeoffs of differential GPS 
further. 
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